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Evidence for an anomalous like-sign dimuon charge asymmetry 
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We measure the charge asymmetry A of hke-sign dimuon events in 6.1 fb~^ of pp collisions recorded 
with the DO detector at a center-of-mass energy y^ = 1.96 TeV at the Fermilab Tevatron collider. 
From A, we extract the like-sign dimuon charge asymmetry in semileptonic 6-hadron decays: A'l = 
-0.00957 ± 0.00251 (stat) ± 0.00146 (syst). This result differs by 3.2 standard deviations from the 
standard model prediction ylgi(5'M) — (— 2.3^Qg) x 10~* and provides first evidence of anomalous 
CP-violation in the mixing of neutral B mesons. 

PACS numbers: 13.25.Hw; 14.40. Nd 



I. INTRODUCTION 



Studies of particle production and decay under the re- 
versal of discrete symmetries (charge, parity and time 
reversal) have yielded considerable insight on the struc- 
ture of the theories that describe high energy phenomena. 
Of particular interest is the observation of CP violation. 
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a phenomenon well established in the K^ and B^ sys- 
tems, but not yet observed for the B^ system, where 
all CP violation effects are expected to be small in the 
standard model (SM) JlJ (See 13] and references therein 
for a review of the experimental results and of the theo- 
retical framework for describing CP violation in neutral 
mesons decays). The violation of CP symmetry is a nec- 
essary condition for baryogenesis, the process thought to 
be responsible for the matter-antimatter asymmetry of 
the universe [J]. However, the observed CP violation in 
the K^ and B^ systems, consistent with the standard 
model expectation, is not sufficient to explain this asym- 
metry, suggesting the presence of additional sources of 
CP violation, beyond the standard model. 

The DO experiment at the Fermilab Tevatron proton- 
antiproton {pp) collider, operating at a center-of-mass 



energy ^/s = 1.96 TeV, is in a unique position to study 
possible effects of CP violation, in particular through the 
study of charge asymmetries in generic final states, given 
that the initial state is CP-symmetric. The high center- 
of-mass energy provides access to mass states beyond the 
reach of the B-factories. The periodic reversal of the DO 
solenoid and toroid polarities results in a cancellation 
at the first order of most detector-related asymmetries. 
In this paper we present a measurement of the like-sign 
dimuon charge asymmetry A, defined as 
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where iV++ and N represent, respectively, the number 
of events in which the two muons of highest transverse 
momentum satisfying the kinematic selections have the 
same positive or negative charge. After removing the con- 
tributions from backgrounds and from residual detector 
effects, we observe a net asymmetry that is significantly 
different from zero. 

We interpret this result assuming that the only source 
of this asymmetry is the mixing of neutral B mesons that 
decay semileptonically, and obtain a measurement of the 
asymmetry A^j defined as 
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where ./V^^ and N^^ 



represent the number of events 
containing two b hadrons decaying semileptonically and 
producing two positive or two negative muons, respec- 
tively. As shown in Appendix |^ each neutral B^ meson 
{q = d, s) contributes a term to this asymmetry given by: 






(3) 



where (fiq is the CP-violating phase, and AMq and Ar,j 
are the mass and width differences between the eigen- 
states of the mass matrices of the neutral B^ mesons. 
The SM predicts the values 0^ = 0.0042 ± 0.0014 and 
(j,^ = _o.096+°;°38 Q. These values set the scale for 
the expected asymmetries in the semileptonic decays of 
B^ mesons that are negligible compared to the present 
experimental sensitivity [l|. In the standard model A^^ 



At,iSM) = (-2.3+-^) X 10 



(4) 



where the uncertainty is mainly due to experimental mea- 
surement of the fraction of B^ mesons produced in pp 
collisions at the Tevatron, and of the parameters control- 
ling the mixing of neutral B mesons. The Pj] semilep- 
tonic charge asymmetry, which constrains the phase (j)^, 
has been measured at e+e~ colliders [21, and the most 
precise results reported bv the BaBar and Belle Collab- 
orations, given in Refs. [J, |5J, are in agreement with the 
SM prediction. Extensions of the SM could produce ad- 
ditional contributions to the Feynman box diagrams re- 
sponsible for B^ mixing and other corrections that can 



provide larger values of (j)q @-0|. Measurements of A^^ 
or (j^q that differ significantly from the SM expectations 
would indicate the presence of new physics. 

The asymmetry A^j is also equal to the charge asym- 
metry Cgj of semileptonic decays of b hadrons to muons 
of "wrong charge" (i.e. a muon charge opposite to the 
charge of the original b quark) induced through B^B^ 
oscillations |1G |: 



<,^ 



TjB -^ 11+ X) - T{B -^ fi-X) 
T{B -^ ii+X) + T{B -^ n-X) 



- A'' 



(5) 



We extract A^j from two observables. The first is the 
like-sign dimuon charge asymmetry A of Eq. ([T]), and the 
second observable is the inclusive muon charge asymme- 
try a defined as 



(6) 



where n+ and ji^ correspond to the number of detected 
positive and negative muons, respectively. 

At the Fermilab Tevatron collider, b quarks are pro- 
duced mainly in bb pairs. The signal for the asymme- 
try A is composed of like-sign dimuon events, with one 
muon arising from direct semileptonic 6-hadron decay 
b — > fi~X [ll|, and the other muon resulting from B'^B^ 
oscillation, followed by the direct semileptonic 5^ me- 
son decay i?° — ?► 5° — ?► ji^X. Consequently the second 
muon has the "wrong sign" due to B^B^ mixing. For 
the asymmetry a, the signal comes from mixing, followed 
by the semileptonic decay i?° — >■ B'^ — ^ fi^ X. The main 
backgrounds for these measurements arise from events 
with at least one muon from kaon or pion decay, or from 
the sequential decay of b quarks fe ^- c — )■ /i+X. For the 
asymmetry a, there is an additional background from di- 
rect production of c-quarks followed by their semileptonic 
decays. 

The data used in this analysis were recorded with the 
DO detector |12h14| at the Fermilab Tevatron proton- 
antiproton collider between April 2002 and June 2009 
and correspond to an integrated luminosity of 6.1 ± 
0.4 fb^^. The result presented in this Article supersedes 
our previous measurement [15J based on the initial data 
set corresponding to 1 fb~^ of integrated luminosity. In 
addition to the larger data set, the main difference be- 
tween these two analyses is that almost all quantities in 
the present measurement are obtained directly from data, 
with minimal input from simulation. To avoid any bias, 
the central value of the asymmetry was extracted from 
the full data set only after all other aspects of the analysis 
and all systematic uncertainties had been finalized. 

The outline of the paper is as follows. In Sec. [TTl we 
present the strategy of the measurement. The detec- 
tor and data selections are discussed in Sec. Illli and in 
Sec. IIVI we describe the Monte Carlo simulations used 
in this analysis. Sections IVllXIIII provide further details. 
Section IXIVI presents the results, Sec. IXVI describes con- 
sistency checks, Sec. IXVII compares the obtained result 



with other existing measurements, and, finahy, Sec lXVlTl 
gives the conclusions. Appendices [K[\E\ provide addi- 
tional technical details on aspects of the analysis. 



II. MEASUREMENT METHOD 

We measure the dimuon charge asymmetry A defined 
in Eq. ([T]) and the inclusive muon charge asymmetry a of 
Eq. ([S]), starting from a dimuon data sample and an in- 
clusive muon sample respectively. Background processes 
and detector asymmetries contribute to these asymme- 
tries. These contributions are measured directly in data 
and used to correct the asymmetries. After applying 
these corrections, the only expected source of residual 
asymmetry in both the inclusive muon and dimuon sam- 
ples is from the asymmetry A^^. Simulations are used to 
relate the residual asymmetries to the asymmetry A^j, 
and to obtain two independent measurements of A^j. 
These measurements are combined to take advantage of 
the correlated contributions from backgrounds, and to re- 
duce the total uncertainties in the determination of A^^. 

The source of the asymmetry a has its nominal ori- 
gin in the semileptonic charge asymmetry of neutral B 
mesons, defined in Eq. ^. However, various detector 
and material- related processes also contribute to n^ . We 
classify all muons into two categories according to their 
origin. The first category, "short" , denoted in the follow- 
ing as "S"' , includes muons from weak decays of b and 
c quarks and r leptons, and from electromagnetic de- 
cays of the short-lived mesons (0, w, 77, p°). The muons 
in the second, "long" , category denoted as "L" , come 
from decays of charged kaons and pious and from other 
processes: charged kaons, pious, and protons not fully 
absorbed by the calorimeter and reaching the muon de- 
tectors ( "punch-through" ) , and false matches of central 
tracks produced by kaons, pious or protons to a track 
segment reconstructed in the muon detector. Thus, the 
L sample contains only the contribution from long-lived 
particles. The total number of muons in the inclusive 
muon sample is 



n = n + n = ns + n^ , 



(7) 



where ns is the number of S muons, and n^ is the number 
of L muons. 

The initial number of observed /i+ (upper signs) or /i^ 
(lower signs) is 



n± ex fsil±as){l±S) + fKil±aK) 
+Uil ± a^) + /p(l ± ap). 



(8) 



In this expression, the quantity S is the charge asymme- 
try related to muon detection and identification, fn is 
the fraction of muons from charged kaon decay, punch- 
through, or false association with a kaon track, and ax 
is their charge asymmetry. This asymmetry is measured 
directly in data as described in Sec. IXIl and therefore, 
by definition, includes the contribution from S. The 



analogous quantities /^r and fp represent the fraction 
of muons from charged pion decay, punch-through or 
false muon association with a pion track, and proton 
punch-through or false muon association with a proton 
track, respectively, while a^ and ap represent the cor- 
responding charge asymmetries. The fraction fp also 
includes a contribution from the association of falsely 
identified tracks with muons. The quantity /bkg = 
nL/{ns+nL) — jK + fix + fp is the L background fraction, 
fs = ns/{ns + nL) = 1 — /bkg is the fraction of S muons, 
and as is related to the semileptonic charge asymmetry 
v4g[, as discussed in Sec. IXIIIl The charge asymmetry a 
can be expressed in terms of these quantities as 



fsias + S) + JKO-K + AOtt + /pflp, 



(9) 



where, because of the small values of 5 and as, only terms 
that depend linearly on the asymmetries are considered. 

The most important background term is /kcik, which 
measures the contribution from kaon decay and punch- 
through. The asymmetry a/f reflects the fact that the in- 
elastic interaction length of the K'^ meson is greater than 
that of the K~ meson ^]. This difference arises from ad- 
ditional hyperon production channels in _ftr~-nucleon re- 
actions, which are absent for their i^"''-nucleon analogs. 
Since the interaction probability of K~^ mesons is smaller, 
they travel further than K~ in the detector material, 
and have a greater chance of decaying to muons, and a 
larger probability to punch-through the absorber mate- 
rial thereby mimicking a muon signal. As a result, the 
asymmetry a^ is positive. Since all other asymmetries 
are at least a factor of ten smaller than ax , neglecting the 
quadratic terms in Eq. Q produces an impact of < 1% 
on the final result. 

In analogy with Eq. ([7]), the number of like-sign 
dimuon events can be written as 



N ^ N++ + N- 



N. 



ss 
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SL 
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LL, 



(10) 



where Nss (Nll) is the number of like-sign dimuon 
events with two S (L) muons, and, similarly, Nsl is the 
number of events with one S and one L muon. A par- 
ticle producing an L muon can be a kaon, pion or pro- 
ton, and, correspondingly, we define the numbers Ngj^ 
with X — K,TT and p. In a similar way, we define N^^ 
with x,y = K,TT,p. The corresponding fractions, de- 
fined per like-sign dimuon event, are f|^ = Ng^/N 
and Fll = Nll/N. We also define Fss = Nss/N, 
FsL = Nsl/N, and Fll = Nll/N. 

The number of observed like-sign dimuon events ^'^ ^.'^ 
(upper signs) or fi~ f.L~ (lower signs) is 

Ar±± ex Fss(l±^5)(l±A)' 

+ Y. ^sl(1±^.)(1±«s)(1±A) 

+ E E F^^lil±A,)il±Ay). (11) 

x,y—K,iT.p; y>x 

The charge asymmetry of Nss events contains the con- 
tribution from the expected asymmetry ^5 that we want 



to measure, and the charge asymmetry A related to the 
detection and identification of muons. The asymmetry of 
the NsL events contains the contribution of background 
asymmetries A^ {x — K,Tr,p) for one muon, and the 
asymmetry (1 ± as){l ± A) for the other muon. The 
asymmetry of Nll events contains the contribution from 
background asymmetries A^ for both muons. By defini- 
tion, the detection asymmetry A is included in the values 
of Ak , Att , and Ap . 

Keeping only the terms linear in asymmetries, the un- 
corrected dimuon charge asymmetry defined in Eq. ([!]) 
can be expressed as 



A = 



N++ - N— 

N++ + N— 

= FssAs + FsLas 

+ (2-Fbkg)A + i^KAK 



F^A^ + FpAp, (12) 



where Fr = F^^ + Ff/ + F^l -h 2F^^ is the total 
number of muons from charged kaon decay or punch- 
through per like-sign dimuon event, and the quantities 
F^ and Fp are defined similarly for charged pious and 
protons. The background fraction -Fbktr is 



-Fbkg — Fk + F.,, + Fp — FsL + 2Fii. 



From Eqs. ([TO]) and ([T3]), it follows that 



FsS+Fy,^„~FLL^l. 



(13) 



(14) 



TABLE I: Fractions of muon candidates in the inclusive muon 
(/^) and in the like-sign dimuon (F^, with two entries per 
event) samples. 



Bin 



Muon pt range (GeV) 






pi- 






1.5 - 2.5 


0.0055 


0.0442 


1 


2.5 - 4.2 


0.1636 


0.2734 


2 


4.2 - 7.0 


0.6587 


0.5017 


3 


7.0 - 10.0 


0.1175 


0.1238 


4 


10.0 - 25.0 


0.0547 


0.0569 



Similarly, the term FkAk is computed as 

4 



FKAK = Y.FlFkc^K, 



(16) 



i=0 



where F* is the fraction of muons in a given pT interval in 
the like-sign dimuon sample. Since the kaon asymmetry 
is determined by the properties of the particle and not 
those of the event, we use the same asymmetry a|c for a 
given Pt interval in both the inclusive muon and the like- 
sign dimuon sample. We verify in Sec. IXVI that the final 
result does not depend significantly on muon -q, nor upon 
kinematic properties of events, luminosity or the mass of 
the /^/J, system. The definition of the muon pT intervals 
and the values of /* and F' are given in Table H] The 
same procedure is applied to all other terms in Eqs. © 
and ((I2]), e.g.. 



As in Eq. ([5]), the largest background contribution in 
Eq. P^ is from the term F^Ak, and all other terms are 
found to be at least a factor of ten smaller. The esti- 
mated contribution from the neglected quadratic terms 
in Eq. p2)) is w 2 x 10~^, which corresponds to sa 4% of 
the statistical uncertainty on A. 

In the following sections, we determine from data all 
the parameters in Eqs. (jH]) and (|12p used to relate the 
measured uncorrected asymmetries a and A to the asym- 
metries as and As- The detection charge asymmetry A 
can differ from S due to differences in the muon trans- 
verse momentum pT and pseudorapidity rj {16.] distribu- 
tions of the like-sign dimuon and inclusive muon data 
samples. For the same reason, we expect the fractions f^ 
in Eq. ^ and F^ in Eq. p^ for x = K,n and p to differ. 
On physics grounds we expect the asymmetries Qx and 
A^ to be identical for any particle of given pT and rj. 

All measurements are performed as a function of the 
muon Pt measured in the central tracker. The range 
oi Pt values between 1.5 and 25 GeV is divided into five 
bins, as shown in TableU] The term Jk^k is obtained by 
the weighted average of the measured values of f}^a}^, i = 
0, 1,2,3,4, with weights given by the fraction of muons 
in a given pT interval, /* , in the inclusive muon sample: 



(2-Fi 



bkg 



fxax = > , fJKa 







I,lJK"'K- 



(15) 



i=0 



F^^5, 



(17) 



As in the case of as, the source of the asymmetry ^5 is 
the charge asymmetry in semileptonic _B-meson decays. 
Thus, two independent measurements of ^^j can be per- 
formed using the inclusive muon and like-sign dimuon 
data samples. The asymmetry as is dominated by detec- 
tor effects, mostly due to the asymmetry arising from the 
different interaction lengths of charged kaons. However, 
As is far more sensitive to the asymmetry ^^j because 
of the definition of A in Eq. ([T]), which has the num- 
ber of like-sign dimuon events, rather than all dimuon 
events in the denominator. Although a weighted average 
of these A\-y measurements can be made, we take advan- 
tage of correlations among backgrounds and asymmetries 
to further improve the precision of A^j through a linear 
combination of A and a. In this combination, which is 
discussed in Sec. IXIVl the detector effects and related 
systematic uncertainties cancel to a large degree, result- 
ing in an improved measurement of ^^j. 



III. DETECTOR AND DATA SELECTION 



The DO detector is described in Refs. [12| - |14j . It con- 
sists of a magnetic central-tracking system that comprises 
a silicon microstrip tracker (SMT) and a central fiber 



tracker (CFT), both located within a 1.9 T supercon- 
ducting solenoidal magnet ^. The SMT has « 800,000 
individual strips, with a typical pitch of 50 — 80 ^m, and a 
design optimized for tracking and vertexing for \ri\ < 2.5. 
The system has a six-barrel longitudinal structure, each 
with a set of four layers arranged axially around the beam 
pipe, and interspersed with 16 radial disks. In the spring 
of 2006, a "Layer 0" barrel detector with 12288 addi- 
tional strips was installed [1J| , and two radial disks were 
removed. The sensors of Layer are located at a ra- 
dius of 17 mm from the colliding beams. The CFT has 
eight thin coaxial barrels, each supporting two doublets 
of overlapping scintillating fibers of 0.835 mm diameter, 
one doublet parallel to the collision axis, and the other 
alternating by ±3° relative to the axis. Light signals are 
transferred via clear fibers to visual light photon counters 
(VLPCs) that have « 80% quantum efficiency. 

The muon system 12] is located beyond the liquid 
Argon-Uranium calorimeters that surround the central 
tracking system, and consists of a layer A of tracking de- 
tectors and scintillation trigger counters before 1.8 T iron 
toroids, followed by two similar layers B and C after the 
toroids. Tracking for JTyj < 1 relies on 10-cm wide drift 
tubes, while 1-cm minidrift tubes are used for 1 < jryj < 2. 

The trigger and data acquisition systems are designed 
to handle the high instantaneous luminosities. Based on 
information from tracking, calorimetry, and muon sys- 
tems, the output of the first level of the trigger is used to 
limit the rate for accepted events to < 2 kHz. At the next 
trigger stage, with more refined information, the rate is 
reduced further to < 1 kHz. These first two levels of trig- 
gering rely mainly on hardware and firmware. The third 
and final level of the trigger, with access to full event 
information, uses software algorithms and a computing 
farm, and reduces the output rate to < 200 Hz, which is 
written to tape. 

The single muon and dimuon triggers used in this anal- 
ysis are based on the information provided by the muon 
detectors, combined with the tracks reconstructed by the 
tracking system. The single muon triggers with the low- 
est pt threshold are prescaled at high instantaneous lu- 
minosity, have a higher average pt threshold than the 
dimuon triggers and cover a smaller range of pseudora- 
pidity than the dimuon triggers. 

In this analysis we select events with one or two muons. 
We therefore first apply track selections, and then require 
either one or two muons. 

Track selection: we select tracks with px in the range 
1.5 < Pt < 25 GeV and \ri\ < 2.2. The upper limit on the 
transverse momentum is applied to suppress the contri- 
bution of muons from W and Z boson decays. To ensure 
that the muon candidate can pass through the detector, 
including all three layers of the muon system, we require 
either px > 4.2 GeV or a longitudinal momentum compo- 
nent \pz\ > 6.4 GeV. The selected tracks have to satisfy 
the following quality requirements: at least 2 axial and 1 
stereo hits in the SMT, and at least 3 axial and 3 stereo 
hits in the CFT. The primary interaction vertex closest 



to this track must contain at least five charged particles. 
This vertex is determined for each event using all re- 
constructed tracks. The average position of the collision 
point in the plane transverse to the beam is measured for 
each run and is used as a constraint. The precision of the 
primary vertex reconstruction for each event is on aver- 
age ss 20 fiia in the transverse plane and ~ 40 fim along 
the beam direction. The transverse impact parameter of 
the selected track relative to the closest primary vertex 
must be < 0.3 cm, with the longitudinal distance from 
the point of closest approach to this vertex < 0.5 cm. 

Single muon selection: the selected track must have a 
matching track segment reconstructed in the muon sys- 
tem, with at least two hits in the layer A chambers, at 
least two hits in the layer B or C chambers, and at least 
one scintillator hit associated with the track. The x^ for 
the difference between the track parameters measured in 
the central tracker and in the muon system must be less 
than 40 (with 5 d.o.f.); the measured time in at least one 
of the scintillators associated with the muon candidate 
must be within 5 ns of the expected time. The muon 
is assigned the charge of the track reconstructed in the 
central tracker. For muon pr < 25 GeV, the fraction 
of muons with mismeasured charge and their contribu- 
tion to the asymmetries are found to be negligible. The 
scintillator timing and the track impact parameter re- 
quirements reduce the background from cosmic rays and 
from beam halo to a negligible level. 

Dimuon selection: The two highest transverse momen- 
tum muons in the event must pass all the selections de- 
scribed above, and be associated to the same interaction 
vertex, applying the same requirements on the transverse 
impact parameter and on the distance of closest approach 
to the primary vertex along the beam axis used in the sin- 
gle muon selection. To remove events in which the two 
muons originate from the decay of the same b hadron, 
we require that the invariant mass of the two muons be 
> 2.8 GeV. 

These requirements define the reference selections, 
which are changed while performing consistency checks 
of the analysis. Unless stated otherwise all figures, tables 
and results in this article refer to these reference selec- 
tions. 

This analysis uses two data samples. The inclusive 
muon sample contains all events with at least one muon 
candidate passing the muon selection and at least one 
single muon trigger. If an event contains more than one 
muon, each muon is included in the inclusive muon sam- 
ple. Such events constitute about 0.5% of the total in- 
clusive muon sample. The like-sign dimuon sample con- 
tains all events with at least two muon candidates of 
the same charge that pass the reference dimuon selection 
and at least one dimuon trigger. If more than two muons 
pass the single muon selection, the two muons with the 
highest Pt are selected for inclusion in the dimuon sam- 
ple. Such events comprise « 0.7% of the total like-sign 
dimuon sample. 

The polarities of the toroidal and solenoidal magnetic 



TABLE II: Weights assigned to the events with different 
solenoid and toroid polarities in the inclusive muon and like- 
sign dimuon samples. 



Solenoid 


Toroid 


Weight 




Weight 


polarity 


polarity 


inclusive muon 


like 


-sign dimuon 


-1 


-1 


0.895 




0.879 


-1 


+1 


1.000 




1.000 


+1 


-1 


0.954 




0.961 


+1 


+1 


0.939 




0.955 



r 0.8 

> 

a 0.6 

a 

:5 0.4 

5 0.2 
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FIG. 1: The normalized muon pr distribution. The points 
correspond to the like-sign dimuon sample and the histogram 
correspond to the inclusive muon sample. The distribution for 
the like-sign dimuon sample contains two entries per event. 



fields are reversed on average every two weeks so that the 
four solenoid-toroid polarity combinations are exposed to 
approximately the same integrated luminosity. This al- 
lows for a cancellation of first order effects related with 
the instrumental asymmetry [15| . To ensure such can- 
cellation, the events are weighted according to the in- 
tegrated luminosity for each dataset corresponding to a 
different configuration of the magnets' polarities. These 
weights are given in Table [TTl 

The normalized pt distributions of muons in the se- 
lected data samples are shown in Fig. [T] Differences 
in these distributions are caused by the trigger require- 
ments. 



IV. MONTE CARLO SIMULATION 

Since almost all quantities are extracted from data, 
the MC simulations are used in only a limited way. The 
simulations of QCD processes used in this analysis are: 

• Inclusive pp collisions containing a minimum trans- 
verse energy i5™'" > 10 or 20 GeV at the generator 
level. 

• Inclusive pp — >■ bbX and pp — >■ ccX final states 
containing a muon, with an additional requirement 



that the b or c quark has transverse momentum 
Pt > 3 GeV, and that the produced muon has px > 
1.5 GeV and \t]\ < 2.1. 

The samples with different ii^™™ are used to study the 
impact of the kinematics of generated events on the pa- 
rameters extracted from the simulation. 

In all cases we use the pythia v6.409 [ij] event gen- 
erator, interfaced to the evtgen decay package [18| and 
the CTEQ6L1 Q^ parton distribution functions. The 
generated events are propagated through the DO detec- 
tor using a GEANT [20] based program with full detector 
simulation. The response in the detector is digitized, and 
the effects of multiple interactions at high luminosity are 
modeled by overlaying hits from randomly triggered pp 
collisions on the digitized hits from MC. The complete 
events are reconstructed with the same program as used 
for data, and, finally, analyzed using the same selection 
criteria described above for data. 



V. MEASUREMENT OF fx, Fr- 

A kaon, pion, or proton can be misidentified as a 
muon and thus contribute to the inclusive muon and the 
like-sign dimuon samples. This can happen because of 
pion and kaon decays in flight, punch- through, or muon 
misidentification. We do not distinguish these individ- 
ual processes, but rather measure the total fraction of 
such particles using data. In the following, the notation 
K ^ II stands for the phrase "kaon misidentified as a 
muon," and the notations n ^ fi and p ~> fi have 
corresponding meanings for pions and protons. In this 
Section we discuss the measurement of fi^ and F^. The 
measurement of the corresponding factors for pions and 
protons and of the asymmetries are discussed in the fol- 
lowing Sections. 

The fraction fn in the inclusive muon sample is mea- 
sured using K**^ -^ K^t:^ decays [lH with K -^ ^x. The 
fraction //f .o of these decays is related to the fraction Jk 
by 



/k-o = eo.fKB.iK*'>), 



(18) 



where R{K*'^) is the fraction of all kaons that result from 
K*^ — > K^Ti^ decays, and Eq is the efficiency to recon- 
struct the pion from the K*^ — )■ K^tt~ decay, provided 
that the K ^ ^j, track is reconstructed. 

We also select Ks mesons and reconstruct K*^ ^■ 
KsTT'^ decays. The number of these decays is 



N{K 



*+ 



Ksn+) = ecN{Ks)RiK*+), (19) 



where R{K*'^) is the fraction of Ks that result from 
K*^ —> KgiT^ decays, and e^ is the efficiency to recon- 
struct the additional pion in the K*^ — ?► Kstt^ decay, 
provided that the Ks meson is reconstructed. We use 
isospin invariance to set 



R{K*'') = R{K*+) 



(20) 



This relation is also confirmed by data as discussed in 
Sec. IVIIII We apply the same kinematic selection cri- 
teria to the charged kaon and Ks candidates, and use 
exactly the same criteria to select an additional pion and 
reconstruct the K*^ — >■ K^it~ and A'*+ — >■ A's7r+ de- 
cays. Therefore we set 



£o — £c 



(21) 



This relation is confirmed by simulation. We assign a sys- 
tematic uncertainty related to this relation, as discussed 
in Sec. Ivml From Eqs. (IIH1)-(I2I1), we obtain 



/ 



K 



N{Ks) 



N{K* 



Ks7r+) 



Ik'O- 



(22) 



We use a similar relation to obtain the quantity Fk of 
K ^> H tracks in the like-sign dimuon sample: 



Fk = 



N{K, 



s) 



N{K*+ ^ Ks7T+) 



Fk'O, 



(23) 



where Fx-o is the fraction of K*'^ -^ K^n^ decays with 
K ^ II in the like-sign dimuon sample. The numbers 
N{Ks) and N{K*+ -> Ksn^) are obtained from the 
inclusive muon sample. 

Since the kaon track parameters must be known to 
reconstruct the K*'^ meson, these measurements of /k 
and Fk require the kaons to decay after being recon- 
structed in the central tracking system. A small num- 
ber of kaon decays occur close to the interaction point, 
so that the muon track is reconstructed by the tracker. 
These muons are counted in the inclusive muon and the 
like-sign dimuon samples, but do not contribute to the 
measurement of the K ^ jj, fraction, because their pa- 
rameters differ significantly from the parameters of the 
original kaon, and they do not produce a narrow K*^ me- 
son peak. The fractions Fk and Jk measured in exclu- 
sive decays are therefore divided by a factor C that cor- 
responds to the fraction of correctly reconstructed kaons 
among aW K ^ ji tracks. This factor is calculated from 
simulation as 



C = 0.938 ±0.006. 



(24) 



Since the mean decay length of kaons in the laboratory 
frame is much longer than the size of the DO detector, 
the value of C is determined mainly by the detector ge- 
ometry, and its value is similar for both K ^ ^i and 
IT —^ fi tracks. Therefore, we use the same coefficient C 
for the computation of the fraction of tt — > /x described in 
Sec. IVIII The difference in this coefficient for kaon and 
pion tracks observed in simulation is taken as the un- 
certainty on its value. The uncertainties from the event 
generation and reconstruction produce a smaller impact 
on this coefficient. 

Details of Kg -^ tt+tt", K*^ -> AT+Tr", and K*+ -> 
KsTT~^ selections and the fitting procedure to measure 
the number of these decays are given in Appendix IB] All 
quantities in Eqs. (|22l) and (1^51) are obtained as a function 



TABLE III: Fractions /k and Fk for different muon pr bins. 
The correspondence between the bin number and the pT range 
is given in Table [H The last line shows the weighted average 
of these quantities obtained with weights given by the fraction 
of muons in a given pr interval /^ (F^) in the inclusive muon 
(dimuon) sample. Only the statistical uncertainties are given. 



Bin 



Ik X 10' 



Fk X 10' 



14.45 ± 1.02 
14.14 ± 0.26 
15.78 ± 0.20 
15.63 ± 0.35 
15.26 ± 0.56 



18.13 ± 4.62 
14.00 ± 1.14 

16.14 ± 0.77 
11.97 ± 1.60 
21.47 ± 2.31 



All 



15.46 ± 0.14 



15.38 ± 0.57 




15 20 25 

Pt(K) [GeV] 



FIG. 2: The fraction of K ^ fi tracks in the inclusive muon 
sample (a) and the like-sign dimuon sample (b), both as a 
function of the pr of the kaon. 



of the measured transverse momentum of the kaon. The 
measured number of K*'^ — >■ A'+tt^ decays with K -^ fi 
in a given px range is normalized by the total number of 
muons in that interval. The fraction Fk'O includes a mul- 
tiplicative factor of two, because there are two muons in 
a like-sign dimuon event, and by definition it is normal- 
ized to the number of like-sign dimuon events. Figure [2] 
and Table IIIII give the resulting fractions /k and Fk for 
different px bins. Only statistical uncertainties are given; 
systematic uncertainties are discussed in Sec. IVIIII 
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VI. MEASUREMENT OF P{-k -^ fJ.)/P{K -> /i) 
AND P(p -^ ^i)/P{K -^ /i) 



The probability P{K — > /i) for a kaon to be misiden- 
tified as a muon is measured using cf) — > K^K~ decays. 
Similarly, we use Kg — )■ 7r"'"7r~ and A — > p7T~ decays |ll| 
to measure the probabilities P{n — > fi) and P{p — >■ ^), 
respectively. In all cases we measure the number N^ 
of decays in which the candidate particle satisfies the 
muon selection criteria defined in Sec. Illli and the num- 
ber of decays iVtr in which the tested particle satisfies the 
track selection criteria. When both kaons (pious) from 
(j> —>■ K^K~^ {Ks —T' TT+TT") satisfy the selection criteria, 
they contribute twice. The details of the event selections 
and of the fitting procedure used to extract the number 
of (j), Ks, and A decays are given in Appendix|B] The ra- 
tio oiNf^{(l)) to Ntr{(p) defines P{K -^ ^)e{^), where e{fi) 
is the efficiency of muon identification. In the same way, 
the ratio of Nf,{Ks) to Ntr{Ks) yields P{tt -^ n)e{fi), 
and the ratio of iV^(A) to A^tr(A) gives the quantity 



P{p -^ ^i)e{^J). The ratios P{tt -^ tJ-)/PiK 
P{p -^ fi)/P{K -^ /i) are obtained from 



P(^^m) _ N^{Ks)/Nt,iKs) 



fi) and 



PiK^fi) 
P{p ^ m) 



N^{A)/Ntr{A) 



PiK~^fi) N^icj))/NtM 



(25) 



Since the initial selection for this measurement requires 
at least one identified muon, we determine all these quan- 
tities in the sub-sample of single muon triggers that con- 
tain at least one muon not associated with the K ^ fi, 
TT — )> /x , or p — ;> /i transitions. 

We measure all these parameters as a function of the 
original particle's transverse momentum. Figures [3] and |4] 
show the ratio P{tt -^ fj.)/P{K -> ^) and P{p -> 
fj,)/P{K — )> fi) respectively, with the mean values av- 
eraged over px determined to be 



P{7T -^ fJ-)/P{K 

p{p ^ ^l)/p{K 



H) =0.540 ±0.029, 
H) =0.076 ±0.021. 



(26) 
(27) 



The dominant uncertainty in Eqs. ([26]) and (|27p stems 
from the limited statistics of data, and the contribu- 
tion of all other uncertainties is much smaller. The 
probability of a pion to be misidentified as a muon 
is much larger than that of proton because the dom- 
inant contribution to this probability comes from the 
7r~ — >■ ^~P decay. The measured ratios ([26]) and ([27]) 
agree well with the results obtained from MC, where we 
obtain P(7r -^ ^i)/P{K -^ A*)(MC) = 0.530 ± 0.011 and 

p{p -^ ^J)/P{K -^ ^JL){MC) = 0.050 ± 0.003. 
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FIG. 3: The ratio P{tx — 5> jj.)/P{K -> /i) as a function of the 
hadron transverse momentum. The horizontal dashed line 
shows the mean value of this ratio. 
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FIG. 4: The ratio P{p — )> p)/P{K -^ /i) as a function of the 
particle transverse momentum. The horizontal dashed line 
shows the mean value of this ratio. 



VII. MEASUREMENT OF ,U, fp, F^, Fp 

The fraction /^ of tt -^ /i tracks in the inclusive muon 
sample can be expressed as 



J-rr Jj 



K 



P{'K — >■ /i) n^ 
P{K^ii)~^' 



(28) 



where the measurement of the fraction /x is described in 
SeclVj that of the ratio P{-k -^ tJ.)/PiK -^ n) in Scc.lVl] 
and the quantities n^ and uk are the mean multiplicities 
of pious and kaons in pp interactions. In a similar way, 
the fraction fp oi p ^>- fi tracks is determined from 



fp = fK 



P{p ^ /i) 



rif 



P{K -^ n) riK 



(29) 



where Up is the average number of protons produced in 
pp interactions. We include in the fraction fp the contri- 
bution from the number Uf of false tracks, reconstructed 
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FIG. 5: The ratio of multiplicities riTr/nK, rip/riK and nt/riK 
as a function of the transverse momentum obtained from 

PYTHIA. 



TABLE IV: Fractions /^, fp, F^, and Fp for different pr bins. 
The correspondence between the bin number and the momen- 
tum range is given in Table|Il The last hne shows the weighted 
averages obtained with weights given by the fraction of muons 
in a given pr interval /^ in the inclusive muon sample. Only 
the statistical uncertainties are given. 



Bin 


/vr X lO'' 


fp X 10^ 


F^ X 10' 


Fp X 10' 


0-1 

2 

3-4 


35.6 ± 4.9 
24.3 ± 1.6 

21.7 ± 3.4 


0.6 ± 0.4 
0.7 ± 0.2 
0.7 ± 0.7 


32.2 ± 5.1 
22.4 ± 1.8 
19.0 ± 3.4 


0.5 ± 0.4 
0.7 ± 0.2 
0.6 ± 0.6 


All 


25.8 ± 1.4 


0.7 ± 0.2 


24.9 ± 1.5 


0.6 ± 0.2 



from random combinations of hits. The impact of false 
tracks on the final result is found to be small and is taken 
into account in the systematic uncertainty. 

The values of uk, n^^, np, and rif are taken from the 
PYTHIA simulation of inclusive hadronic interactions. We 
count the number of particles satisfying the track selec- 
tion criteria in the simulated interactions, and obtain the 
dependence of the ratios riT^/riK, np/riK and nf/tiK on 
the particle px shown in Fig. [51 

Both Jtt and fp are measured as a function of the par- 
ticle Pt- However, they are poorly defined in the first 
and last bins due to low statistics. Therefore, we com- 
bine these quantities for bins and 1 and for bins 3 and 
4. 

Figure [6] and Table IIVI provide the measured fractions 
/^ and fp for different pr bins. Only statistical uncer- 
tainties are given. The systematic uncertainties related 
to these quantities are discussed in Sec. I Villi 

The fractions -FV and Fp in the like-sign dimuon sample 
are determined in a similar way: 



F^ - Fj 



K 



F„ 



F, 



K 



P{K -^ti)NK' 
P{p -^p.) Np + Nf 
P(K ^ ^) Nk 



(30) 
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FIG. 6: The fraction of (a) tt — >■ /i tracks and (b) p — > 
/i tracks in the inclusive muon sample as a function of the 
track transverse momentum. The horizontal dashed lines 
show the mean values of these fractions. 



where the quantities Nx, -/V^r, Np, and Ni represent 
the average numbers of kaons, pions, protons and false 
tracks for events with two identified muons with the same 
charge. The simulation shows that the ratio Nt^/Nk can 
be approximated as 



— = (0.90 ±0.05) — . 
Nk UK 



(31) 



The main uncertainty in this value is due to the simula- 
tion of pion and kaon multiplicities in pp interactions and 
is discussed in Sec. I Villi The ratio Np/Nx is also consis- 
tent with the factor given in Eq.|3T] The value of N-j^/Nk 
is smaller than that oin-,^lnx because the main contribu- 
tion in the sample with one identified muon comes from 
semileptonic decays of b and c quarks, which usually also 
contain at least one kaon. Since the number of simulated 
events with one identified muon is small, we obtain the 
ratios N.^/Nk and {Np -f N{)/Nk using the approxima- 
tion of Eq. ((3T|) . i.e., multiplying the quantities n-^/nK 
and {np + nf)/nK by the factor 0.90 ±0.05. Figure[7|and 
Table IIVI give the fractions i^^ and Fp for different pt 
bins. Only the statistical uncertainties of the simulation 
are given. The systematic uncertainties related to these 
quantities are discussed in Sec. I Villi As in the case of 
/tt and fp, the mean value of these quantities are used in 
bins and 1 and for bins 3 and 4. 
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FIG. 7: The fraction of (a) tt — >■ /i tracks and (b) p — > 
ji tracks in the hke-sign dimuon sample as a function of 
the track transverse momentum. The horizontal dashed lines 
show the mean values of these fractions. 



SYSTEMATIC UNCERTAINTIES OF 
BACKGROUND FRACTIONS 



VIII. 



We use Eqs. (PH)) and (pi]) to derive the fractions /k 
and Fk, and verify the validity of Eq. (PT|) in simulation, 
and find that 



(sc/eo) 



MC 



0.986 ±0.029, 



(32) 



where the uncertainty reflects only the statistics of the 
Monte Carlo. The ratio R{K*+)/R{K*°) measured in 
simulation is 



R{K*+) / R{K*°) = 0.959 ± 0.035. 



(33) 



The validity of Eq. (|20p in simulation relies mainly on the 
assumptions used in the fragmentation and hadronization 
processes in the event generator. To confirm the valid- 
ity of Eq. (PO)) . we use the existing experimental data on 
K'^ , Ks, K*'^, and K*~^ multiplicities in jets, which were 
obtained at e+e" colliders at different center-of-mass en- 
ergies [2]. From these data we obtain: 

R{K*+)/R{K*°) ^ 1.039 ±0.075. 



(34) 



The simulation and data are consistent, and we assign a 
relative uncertainty of 7.5% to both Jk and Fk from the 
assumption of Eq. pp)) . We also assign an uncertainty of 
4% due to the fitting procedure used to extract the num- 
bers of _ftr*+ and K*'^ events. This uncertainty is obtained 



by varying the background parametrization and the fit- 
ting range. Since the same background model is used to 
obtain the number of K*'^ events, both in the inclusive 
muon and the like-sign dimuon samples, this uncertainty 
is taken to be the same for Jk and Fk- Adding all con- 
tributions in quadrature, and including the uncertainty 
in Eq. (I24|) . we find a relative systematic uncertainty of 
9.0% in Jk and Fk, with a 100% correlation between the 
two. 

We assign an additional uncertainty of 2.0% on Fk 
due to the description of the background in the inclu- 
sive muon and like-sign dimuon events. This uncertainty 
is estimated by varying the background parametrization 
and range used for fitting, and by comparing with the 
results of the alternative fitting method presented in Ap- 
pendix |e1 

We use the ratio N{Ks)/N{K*+ — > Kstt^) to measure 
both fK and Fk, which is equivalent to the statement 
that the ratios Fk*o/Fk and fK'°/ Ik are identical, i.e., 
that the fraction of kaons originating from K*'^ is the 
same in the inclusive muon and in the like-sign dimuon 
samples. This is validated in simulation with an uncer- 
tainty of 3% due to the statistics of the simulation, which 
we assign as an additional systematic uncertainty to the 
fraction Fk- 

The uncertainty on the background fractions /^, i^^, 
fp, and Fp have an additional contribution from the ra- 
tios of multiplicities riT^/nK and Up/riK extracted from 
the simulation. To test the validity of the simulation, we 
measure the multiplicity of kaons in the inclusive muon 
sample. We select events with one reconstructed muon 
and at least one additional charged particle that satisfies 
the track selection criteria. We determine the fraction 
of kaons among these tracks using the same method as 
in Sec. |Vl i.e., we find the fraction of tracks from the 
j^*Q _j, j{+.jY~- decay and convert this into the fraction of 
kaons. We compare the kaon multiplicity in data using 
this method with that measured in the simulation, and 
we find that they agree within 10%. Since part of this 
difference can be attributed to the uncertainties from the 
assumptions of Eqs. ([^0]) and (|2ip . and part is due to the 
fitting procedure described above, we find that the un- 
certainty of the kaon multiplicity in the simulation does 
not exceed 4%, and we assign this uncertainty to both 
quantities ny^/nK and Up/riK- We also assign this 4% 
uncertainty to the Eq. (|3ip used to derive the values of 
NJNk and Np/NK- 

Any falsely reconstructed track identified as a muon is 
treated in the analysis in the same way as a proton. We 
check the impact of this approach by completely remov- 
ing the contribution of false tracks, or by increasing their 
contribution by a factor of ten, and the final value of A'^^ 
changes by less than 0.00016. We include this difference 
as the systematic uncertainty on the contribution from 
false tracks. 
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TABLE V: Fractions fs, /k, f-w, and fp + fi, measured in data 
and in simulation (MC). Only the statistical uncertainties on 
these measurements are shown. 

fs X 10^ /a- X 10^ U X 10'' ifp + /f) X 10^ 



Data 58.1 ± 1.4 
MC 59.0 ± 0.3 



15.5 ± 0.2 
14.5 ± 0.2 



25.9 ± 1.4 
25.7 ± 0.3 



0.7 ± 0.2 
0.8 ± 0.1 



IX. MEASUREMENT OF fs, Fss 

We use the measurements of the fractions of back- 
ground muons in the "long" category, obtained in the 
previous sections to evaluate the fraction of muons in 
the "short" category. In the inclusive muon sample the 
fraction fs is determined as 



/s — 1 - Ik - fit — fp 



(35) 



We check through simulation that the contribution from 
all other sources to the inclusive muon sample, such as 
Kl — >■ TT/.t^' decays or the semileptonic decays of hypcr- 
ons, is negligible. The muons from t — > i-w^j_Vr arc in- 
cluded by definition in the fs- The fraction fs is mea- 
sured separately in each muon pT bin and then a weighted 
average is calculated with weights given by the fraction 
of muons in a given p^ interval fl in the inclusive muon 
samples. From data, we obtain 



fs = 0.581 ± 0.014 (stat) ± 0.039 (syst). 



(36) 



where the systematic uncertainty comes from the uncer- 
tainty on the background fractions described in Sec. I Villi 

To check the procedure for determining the back- 
ground fractions, the composition of the inclusive muon 
sample in data is compared to that from simulation in Ta- 
ble |Vl where only statistical uncertainties for both data 
and simulation are shown. The agreement between data 
and simulation is very good, and the remaining differ- 
ences are within the assigned systematic uncertainties. 
Although the values given in Table |V] for data and for 
simulation are not independent, some, such as fx and 
P{-K — 7^ ^)/P{K — >• fj.) used to derive /tt, are measured 
directly in data. As a consequence, this result can be 
used as an additional confirmation of the validity of our 
method. 

The background fractions Fk, F-^, and Fp are obtained 
from the same weighted average used for Eq. p6p of the 
quantities measured in each px interval i, starting from 
the values given in Sees. IVl and IVIII ff^' is used as weight 
instead of /^). Using Eq. P^ . we obtain 



i^bkg = 0.409 ± 0.019 (stat) ± 0.040 (syst). 



(37) 



To evaluate the fraction Fss in the like-sign dimuon 
sample, we take into account that, in some events, both 
muons belong to the L category. The fraction of these 
events in all events with at least one L muon is measured 
in simulation and found to be 

Fll 



The uncertainty in Eq. (|38|) includes the 4% systematic 
uncertainty related to the multiplicity of different parti- 
cles in the simulation, as discussed in Sec. IVIIII Using 
Eqs. dni), dSZl), and dMl, we obtain 

Fll = 0.074 ± 0.003 (stat) ± 0.008 (syst), (39) 

and, finally, from Eqs. (ITil) . (jSZ]), and (P^l we obtain 

Fss = 0.665 ± 0.016 (stat) ± 0.033 (syst). (40) 

X. MEASUREMENT OF 5 

Table III of Ref. [1^1 gives a complete list of the contri- 
butions to the dimuon charge asymmetry that are caused 
by detector effects. The largest of these effects is « 3%. 
The reversal of magnet polarities is a characteristic of the 
DO experiment that allows the cancellation at first order 
of these detector effects, reducing any charge asymm etry 
introduced by the track reconstruction considerably [15| . 

Higher-order effects result in a small residual recon- 
struction asymmetry at the 10~^ level. This asymmetry 
is measured using J/V' — J- /^^M^ decays. We select events 
that pass at least one dimuon trigger and have at least 
one identified muon and one additional particle of op- 
posite charge that satisfies the track selection criteria of 
Sec, mil We veriiy in Appendices [Q and ID] that the track 
reconstruction and trigger selection do not introduce an 
additional charge asymmetry. The residual asymmetry 
is measured as a function of the muon transverse mo- 
mentum. The probability to identify a muon with charge 
Q = ±1 and pt corresponding to bin i is denoted by 
P,;(l + Q5i)^ where Pi is the mean probability for pos- 
itive and negative muons, and 5i is the muon detection 
asymmetry we want to measure. The probability of iden- 
tifying the second muon with pT in bin j , provided that 
the first muon has pt in bin i, is denoted by P*(l -t- Qdj). 

The number of events Nij with a positive muon in bin 
i and negative muon in bin j is 



N,,^NP,{l + 6,)P;{l-6j), 



(41) 



where N is the total number of selected J/V' — >• /^"""M" 
decays. The number of events with only one selected 
muon of charge Q is 



where the sum extends over the five transverse momen- 
tum intervals. 

The probabilities Pi and P!- are not independent. From 
Eq. (|4ip . we have the following normalization condition 



FsL + Fll 



0.220 ±0.012. 



(38) 



P.P'=P,P' 



J — ^ J^ i 



(43) 
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In addition, since the total probability to identify the 



muon is Ptot = J2i 
tion condition 



_(j Pi , we get the following normaliza- 



(Ptot)' = E p^Pj 



(44) 



Experimentally we measure the quantities Na, E^ , A^j 
{i < j), Si, and A^, which can be expressed as 



whe: 
pi \ 

sum/ 

Since the number of measured quantities is greater 
than the number of unknowns. Pi, P^, and Si can be 
obtained from Eqs. (I43H45|) by minimizing the x^ of the 
difference between the observed and expected quantities. 
The quantities Na, T^ij, Ay (i < j), E.^, and A^ are 
obtained from fits to the J/ip mass peak in the dimuon 
invariant mass distribution M{fjr^ ii~) in each of the five 
Pt bins. The J/i/' signal is described by the sum of two 
Gaussians. An additional Gaussian is included to take 
into account the contribution from the if)' . The back- 
ground is parametrized by a third degree polynomial. 
The mean position and R.M.S. of all Gaussian functions 
in the fit of Ay and A^ are fixed to the values obtained in 
the fit of the corresponding quantities S^ and S^. Exam- 
ples of the fits are shown in Figs. |8] and [9l The values of 
5i obtained as a function of the muon px are given in Ta- 
ble lVII and are shown in Fig.[TOl The correlations between 
values of 5i in different bins are given in Table I VIII The 
weighted average for the residual muon asymmetry in the 
inclusive muon and the like-sign dimuon samples, calcu- 
lated using weights given by the fraction of muons in a 
given Pt interval /* (P') in the inclusive muon (dimuon) 
sample, are given respectively by 



5= Etoi^-^- 



-0.00076 ±0.00028, (46) 

-0.00068 ±0.00023, (47) 



where only the statistical uncertainties are given. The 
correlations among different 5i are taken into account. 
These small values of the residual muon reconstruction 
asymmetry are a direct consequence of the regular rever- 
sal of the magnets polarities in the DO experiment. 
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FIG. 8: The /i^/^ invariant mass distributions used to obtain 
(a) E23 = iV23 + N'j,2 and (b) A23 = iV23 - A^32. The solid 
line presents the result of the fit; the dashed line shows the 
non-resonant background contribution. 
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XI. MEASUREMENT OF aK, a^, fflp 

The largest detector-related charge asymmetry is pro- 
duced by -JC — >■ /I tracks. It is caused, as discussed in 



FIG. 9: The /i /i invariant mass distributions used to obtain 
(a) E2 = ^^2+ + A'^2- and (b) A2 = N2+ - N2-. The solid 
line presents the result of the fit; the dashed line shows the 
non-resonant background contribution. 



15 



TABLE VI: Muon reconstruction asymmetry Si for difTerent 
muon pt bins. The correspondence between the bin number 
and the pT range is given in Table IH Only the statistical 
uncertainties are given. 

Bin Si 



-0.00203 ± 0.00194 
-0.00045 ± 0.00059 
-0.00130 ± 0.00048 
+0.00075 ± 0.00125 
+0.00162 + 0.00230 



TABLE VIL Correlation coefficients among values of Si in 
different bins. The correspondence between the bin number 
and the pr range is given in Table IH 



Bin 





1 


2 


3 


4 





+1.000 


-0.189 


-0.155 


+0.024 


-0.051 


1 


-0.189 


+1.000 


-0.449 


-0.117 


-0.059 


2 


-0.155 


-0.449 


+1.000 


-0.242 


-0.124 


3 


+0.024 


-0.117 


-0.242 


+ 1.000 


-0.006 


4 


-0.051 


-0.059 


-0.124 


-0.006 


+1.000 



Sec. nil by the difference between the K^N and K^N in- 
teraction cross sections [2| , resulting in a positive charge 
asymmetry of muons coming from kaon decay or punch- 
through. 

The asymmetry ax oi K -^ ii tracks is measured 
directly in data using K*° — )• K+tt~ and </> — >■ K+K~ 
decays. In both cases we select candidates with K ^■ 
/i tracks, in the entire inclusive muon sample. We calcu- 
late separate mass distributions for positive and negative 
K —> fi tracks, and fit the sum and the difference of these 
distributions to extract the quantity A^^ , corresponding 
to the difference in the number of K*^ or (p meson de- 
cays with positive and negative K ^f ^i tracks, and the 
quantity T^k, corresponding to their sum. The selection 
of events and the fitting procedure used to extract the 
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FIG. 10: Muon reconstruction asymmetry as a function of the 
muon Pt- 



number of signal decays are described in Appendix |B] 
The asymmetry gk is measured as: 



a-K = C/^k/^k, 



(48) 



where the coefficient C is the fraction of correctly recon- 
structed kaons among all iiT — > /i tracks as in Eq. (|24p . 
In this measurement of qk, we require that the kaon 
decays after having been reconstructed in the tracking 
system, since its track parameters must be measured in 
order to reconstruct the K*^ or meson. However, the 
K ^ ^ tracks in the inclusive sample also include kaons 
decaying before being reconstructed in the tracker. Since 
the kaon asymmetry is caused by the interactions of kaons 
with the material of the detector, and the amount of ma- 
terial near the interaction point is negligible, the kaons 
decaying before being reconstructed by the tracker do 
not produce any significant asymmetry. They contribute 
only in the denominator of Eq. (|48l) . The factor C takes 
into account the contribution of these tracks. Its numeri- 
cal value is given in Eq. (l24t . It should be noted that this 
factor cancels in the products /ko-k, etc., since both fx 
and ax are measured using the correctly reconstructed 
K ^ fi tracks. 

Figure [TlTa) shows the value of a^iK*^) measured 
in K*'^ -^ A'+7r~ decay as a function of the pt of the 
K ^1' fi track. The asymmetry in (/) — >■ K^K'^ decays 
has to be corrected for the charge asymmetry of 



the second kaon track a}^'^'^'^: 



aK{cj)) 



^K 



^ track 



(49) 



The kaon track reconstruction asymmetry a^^"^*^ is dis- 
cussed and measured as a function of kaon momentum 
in [22| using the decay D*^ -^ D^n'^ with D" ^■ 
K~ jX^v^ and is taken from that article. It is convoluted 
with the Pt distribution of the second kaon for each bin 
of the K -^ fJ.pT- Figure fTlT b) shows the resulting asym- 
metry axiip)- 

The two measurements of uk are consistent. The 
X^/d.o.f . for their difference is 5.40/5. Therefore they can 
be combined and the resulting asymmetry gk is shown 
in Fig. [T^l and in Table IVIIII Due to the requirement of 
Pt > 4.2 GeV or \pz\ > 6.3 GeV, the first two bins in 
Fig. [T2] correspond to muons that traverse the forward 
toroids of the DO detector. These muons have a larger 
momentum and a longer path length before the calorime- 
ter than central muons. As a result, ax drops at low pT- 

The asymmetry Ott of tt — >■ /x tracks and the asymme- 
try ftp of p — >■ /x tracks are expected to be much smaller. 
We measure these asymmetries using A'5 — > tt+tt^ and 
A -+ pTT~ decays, respectively. The details of the Ks 
and A selections are given in Appendix [Bj The tech- 
nique used to measure the asymmetry is the same as in 
UK measurement. The same factor C is used to measure 
the asymmetry ajr from Kg — > tt+tt" decays. The un- 
certainty in C takes into account its difference for kaon 
and pion tracks. Since the proton is stable, this fac- 
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TABLE VIII: Asymmetries ax, cltt, and Op for different pr 
bins. Tfie correspondence between the bin number and px 
range is given in Table [H Tlie last line shows the mean asym- 
metries averaged over the inclusive muon sample. Only the 
statistical uncertainties are given. 



Bin 



cs 



aK 



+0.0526 ± 0.0242 



+0.0027 ±0.0021 -0.104 + 0.076 



1 +0.0424 ± 0.0027 

2 +0.0564 ± 0.0013 +0.0013 + 0.0012 +0.028 + 0.035 

3 +0.0620 + 0.0032 



+0.0620 + 0.0048 



+0.0164 + 0.0044 +0.077 + 0.055 



AH +0.0551+0.0011 +0.0025 + 0.0010 +0.023 + 0.028 
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FIG. 11: The asymmetry aK measured with (a) K*^ — >■ 
K^-k' and (b) — >■ K^ K~ decays as a function of the 
K ^ fi pt- 
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FIG. 13: The asymmetry (a) a^ and (b) Up as a function of 
the Pt of the pion and proton, respectively. 



tor is not used in the computation of the asymmetry of 
p — )■ /x tracks. 

The asymmetries a^ and ap are shown in Fig. [13] as 
a function of the px of the tt ^ /i and p — >■ /i tracks, 
respectively. The values of these asymmetries and their 
averages are Usted in Table IVIIII for different px bins. 
We use the mean value of these quantities in bins and 
1 and in bins 3 and 4 since the statistics available in the 
first and last bin are not sufficient to perform separate 
measurements . 

The asymmetries Ak, A.^ and Ap are obtained from 
ttK , flir and Op using Eq. [12] (and analogous relations for 
pions and protons). 
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FIG. 12: The combined asymmetry ax as a function of pr. 



XII. CORRECTIONS DUE TO BACKGROUND 
ASYMMETRIES 



The corrections for the asymmetries of the background, 
obtained from Tables HI IIIIl IIVI and I Villi are summarized 
in Tables IIXI and [X] The values Jko-K: FkAk, etc., 
arc computed by averaging the corresponding quantities 
with weights given by the fraction of muons in a given px 
interval fl{Fl) in the inclusive muon (dimuon) sample, 
see Eqs. (|l5l) and (fT6|) . We use the mean value of Z^, F^, 
/p, Fp, a^, and ap in bins and 1 and in bins 3 and 4 
as the statistics available in the first and last bin are not 
sufficient to perform separate measurements. 



17 



TABLE IX: Corrections due to background asymmetries 
fKO-K, fnin, aud fpUp for different pt bins. The last line 
shows the weighted averages obtained using weights given by 
the fraction of muons in a given pr interval /^ in the inclusive 
muon sample. Only the statistical uncertainties are given. 



Bin 



fKaK X 10^ 



/vrClTr X 10 



fpttp X 10 



+0.760 ± 0.353 

1 +0.600 ± 0.040 

2 +0.889 ± 0.023 

3 +0.968 ± 0.054 

4 +0.946 ± 0.081 



+0.095 ± 0.076 
+0.033 ± 0.030 
+0.337 + 0.109 



-0.061 + 0.060 
+0.020 + 0.026 
+0.053 + 0.067 



AH +0.854 + 0.018 +0.095 + 0.027 +0.012 + 0.022 



TABLE X: Corrections due to background asymmetries 
FkAk, FttAtt and FpAp for different pr bins. The last line 
shows the weighted averages obtained using weights given by 
the fraction of muons in a given px interval F^ in the dimuon 
sample. Only the statistical uncertainties are given. 



Bin 



FkAk X 10' 



F^A.„ X 10' 







+0.953 + 0.501 



+0.086 ± 0.069 



-0.056 + 0.054 



1 +0.594 + 0.061 

2 +0.910 + 0.048 +0.030 + 0.027 +0.019 + 0.024 

3 +0.741 + 0.106 



+1.332 + 0.176 



+0.294 + 0.098 +0.046 + 0.058 



AH +0.828 + 0.035 +0.095 + 0.025 +0.000 + 0.021 



XIII. ASYMMETRIES as AND As 

In the absence of new particles or interactions, the only 
non-instrumental source of the asymmetries as and As is 
the semileptonic charge asymmetry A^j given by Eq. ([SJ . 
Both as and As are proportional to A^j, through the 
coefficients 



Cb = 



asMsh 
As/Al,, 



(50) 
(51) 



which are determined from simulation. 

The decays producing an S muon in the inclusive muon 
sample, and their weights relative to the semileptonic 
decay h — >■ ^X [ll| , are listed in Table IXII All weights 
are computed using simulated events. The main process, 
denoted as Ti, is the direct semileptonic decay of a 6 
quark. It includes the decays b — > ^X and b — > rX, 
with T — >■ ^X. The weights wia and wih for semileptonic 
decays of B mesons with and without oscillations are 
computed using the mean mixing probability 



XO = f'dXdQ + /sXsO, 



(52) 



where /^ and /^ are the fractions of B^ and B^^ mesons 
in a sample of semileptonic i?-meson decays, and Xdo and 
XsO are the B^ and B^ mixing probabilities integrated 
over time. We use the value Xq_= 0.147± 0.011 measured 
at the Tevatron and given in [2, [23 • The second process 
T2 concerns the sequential decay & — >■ c — >■ /iX. For 



TABLE XI: Heavy quark decays contributing to the inclu- 
sive muon and like-sign dimuon samples. Abbreviation "nos" 
stands for "non-oscillating," and "osc" for "oscillating." All 
weights are computed using the MC simulation. 



Process 



Weight 



Ti b^ fi'X 

Tia b ^ II- X (nos) 

Tib b ^f h ^f jjT X (osc) 

Ta b^ c^ /i+X 

T2a b ^ c^ ji^X (nos) 

T2b 5 — > 6 — ^ c — ^ /i+X (osc) 



Wl = 1. 

Wla = (1 — Xo)™l 

Wib = Xowi 

W2 = 0.113 + 0.010 

W2a = (1 - Xo)W2 
W2b = X0'W2 



Ti b^ ccq with c -^ (1+ X or c -> ^'X W3 ^ 0.062 + 0.006 
T4 77^0;, p°, 0(1020), J/V'.V'' ^M^M" W4 = 0.021 + 0.001 
Ts bbcc with c ->• /i+X or c ^ fi' X W5 = 0.013 + 0.002 



Te cc with c — >■ /i X or c — > /i X 



W(i 



0.660 + 0.077 



simplicity we use the same value of xo to compute the 
weights W2a and W2b of non-oscillating and oscillating 
sequential decays 6 — !• c — ?► fiX. The process T3 is the 
decay of a 6 hadron to a cc pair, with either the c or c 
quark producing a muon, while T4 includes the decays of 
short-lived mesons 77, ui, p", (j){1020), J/ip, and tp' to a 
^+/.j^ pair. We take into account both the decays of b 
hadrons to these particles and their prompt production. 
The process T5 represents four-quark production of bbcc 
with either the c or c quark decaying to a muon. The 
decays of 6 or 6 quark to a muon in this process and 
the four-quark production of bbbb are taken into account 
through processes Ti, T2, and T3. Finally, the process 
Tg involves cc production followed by c — ?► fiX decay. 
We separate the processes T5 and Tg because only T5 
contributes to the like-sign dimuon sample, while both 
T5 and Tq contribute to the inclusive muon sample. 

The uncertainty in the weights of different processes 
contains contributions from the uncertainty in the mo- 
mentum of the generated b hadrons and from the un- 
certainties of branching fractions for 6-hadron decays. 
We reweight the simulated 6-hadron momentum to get 
agreement of the of muon momentum spectrum in data 
and in MC, and the difference in weights is assigned as 
the systematic uncertainty on the momentum distribu- 
tion. The uncertainties in the inclusive branching frac- 
tions B — 7> fiX, B — > cX and B — > cX taken from 2] are 
propagated into the uncertainties on the corresponding 
weights. We assign an additional uncertainty of 10% to 
the weights w^ and wq due to the uncertainties on the 
production cross sections of cc and bbcc processes. 

Among all processes listed in Table IXll the process Tu, 
is directly related to the semileptonic charge asymmetry 
ylgj (see Appendix |X] for details) . The process T2b pro- 
duces the flavor-specific charge asymmetry Afs- We set 
Afs — — ^si, where the negative sign appears because the 
charge of the muon in the process T26 is opposite to the 
charge of the muon in the process Tk,. No other process 
contributes to the charge asymmetry and therefore they 

The coefficient Cb is found 



just dilute the value of ^^j 
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from: 



Cb 



Wlb - W2b 



Wi + W2 + W3 + W4 + W5 + Wq 



0.070 ±0.006. 



(53) 

The computation of the coefficient Cb is more comph- 
cated. One of the selections for the hke-sign dimuon sam- 
ple requires that the invariant mass of the two muons be 
greater than 2.8 GeV. This requirement suppresses the 
contribution from processes in which both muons arise 
from the decay of the same quark. The probability that 
the initial b quark produces a /i" is 



Pb ex Wia + W2b + 0.5(W3 + k;4 + W5), 



(54) 



where we apply the coefficient 0.5 because processes T3, 
T4, and T5 produce an equal number of positive and neg- 
ative muons. The probability that the accompanying b 
quark also produces a fi^ is 



Pfc ex Wlb + W2a + 0.5(W3 +W4 + W5). 



(55) 



The total probability of observing like-sign dimuon events 
from decays of heavy quarks is 



Ptot (X PbP-b- 



(56) 



The probability of processes contributing to the charge 
asymmetry of dimuon events is 

Pas OC Wib[wia+0.5{W3 +W4 + IU5)] - 

W2b[w2a+0.5{W3+W4+W5)]. (57) 

The coefficient Cb is obtained from the ratio 

Cb =- Pas/Ptot = 0.486 ± 0.032. (58) 

This relation assumes that the processes producing the 
two muons are independent and is verified by calculating 
the coefficient Cb in simulated like-sign dimuon events. 
We exclude the process Tg with cc pair production, since 
the mixing probability of D^ meson is small and these 
events do not contribute significantly to the like-sign 
dimuon sample. We count the number of direct-direct 
6-quark decays, Ndd, of direct-sequential decays, Nds, 
of sequential-sequential decays, Ngs, of direct-random 
events, Ndr ("random" includes processes T3, T4, and 
T5), of sequential-random decays, Nsr, to obtain 



Cb = 



N, 



dd 



Nss + XQiNdr ~ Nsr) 



Nls 



= 0.448±0.071, (59) 



where Nis is the total number of like-sign dimuon events. 
This result agrees well with the value in Eq. (|58p . The 
uncertainty of this method is larger because of the small 
statistics of simulated like-sign dimuon events. 



XIV. ASYMMETRY A^i 

The uncorrected asymmetries a and A are obtained by 
counting the number of events of each charge in the inclu- 
sive muon and the like-sign dimuon samples, respectively. 



TABLE XII: Sources of uncertainty on yl,^ in Eqs. ((62]), (f63l) . 
and (|65p . The first eight rows contain statistical uncertainties, 
the next three rows contain systematic uncertainties. 



Source 


5o{Al,)m 


5a{Al{)m 


5a{Al,)m 


yl or a (stat) 


0.00066 


0.00159 


0.00179 


fK or Fk (stat) 


0.00222 


0.00123 


0.00140 


P{7V ^ ^,)/P{K ^ /i) 


0.00234 


0.00038 


0.00010 


P{p ^ f,)/P{K ^ m) 


0.00301 


0.00044 


0.00011 


Ak 


0.00410 


0.00076 


0.00061 


A^ 


0.00699 


0.00086 


0.00035 


Ap 


0.00478 


0.00054 


0.00001 


S or A 


0.00405 


0.00105 


0.00077 


fK or Fk (syst) 


0.02137 


0.00300 


0.00128 


IT, K, p multiplicity 


0.00098 


0.00025 


0.00018 


Cb or Cb 


0.00080 


0.00046 


0.00068 


Total statistical 


0.01118 


0.00266 


0.00251 


Total systematic 


0.02140 


0.00305 


0.00146 


Total 


0.02415 


0.00405 


0.00290 



In total, there are 1.495 x 10^ muons in the inclusive muon 
sample, and 3.731 x 10^ events in the like-sign dimuon 
sample. We obtain 



a = +0.00955 ±0.00003, 
A = +0.00564 ±0.00053. 



(60) 
(61) 



The results obtained in Sees. IVHXIIII are used to cal- 
culate the asymmetries as and As from these values, 
which are then used to evaluate the charge asymmetry 
for semileptonic B meson decays. 

The asymmetry A^j, extracted from the asymmetry a 
of the inclusive muon sample using Eqs. ([9|) and ([50]) . is 



Al^ = +0.0094 ± 0.0112 (stat) ± 0.0214 (syst). 



(62) 



The contributions to the uncertainty on this value are 
given in Table IXIII Figure [T4l(a) shows a compari- 
son of the asymmetry a and the background asymmetry 
flbkg = fs^ ± IkIk + fna-K + fpO-p, as a function of the 
muonpT- There is excellent agreement between these two 
quantities, with the x^/d.o.f. for their difference being 
2.4/5. Figure [Mfb) shows the value of fs^s = a — abkg, 
which is consistent with zero. The values a and flbkg are 
given in Table IXlTTl This result agrees with expectations, 
since the value of the asymmetry a should be determined 
mainly by the background, and the contribution from 
Agj should be strongly suppressed by the small factor of 
C6 = 0.070 ±0.006. 

The consistency of ^^j with zero in Eq. (|62p and the 
good description of the charge asymmetry a for different 
values of the muon px shown in Fig. [T3] constitute an 
important tests of the validity of the background model 
and of the analysis method discussed in this article. 

The second measurement of the asymmetry A^^, ob- 
tained from the uncorrected asymmetry A of the like-sign 
dimuon sample using Eqs. (|12p and (|5ip . is 



A'' 



-0.00736 ± 0.00266 (stat) ± 0.00305 (syst). (63) 
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FIG. 14: (a) The asymmetry Hbkg (points with error bars) as 
expected from our measurements of the fractions and asym- 
metries of the background processes is compared to the mea- 
sured asymmetry a of the inclusive muon sample (shown as 
histogram, since the statistical uncertainties are negligible). 
The asymmetry from CP violation is neghgible compared to 
the background in the inclusive muon sample; (b) the differ- 
ence a — flbkg. The horizontal dashed line shows the mean 
value of this difference. 



TABLE XIII: The measured asymmetry a and the expected 
background asymmetry abkg in the inclusive muon sample for 
different pr bins. For the background asymmetry, the first 
uncertainty is statistical, the second is systematic. 



bin 


a X 10^ 


abkg X 10^ 



1 
2 
3 

4 


0.324 ± 0.036 
0.582 ± 0.007 
0.978 ± 0.003 
1.193 ± 0.008 
1.339 ± 0.011 


0.693 ± 0.379 ± 0.632 
0.611 ± 0.109 ±0.072 
0.865 ± 0.054 ± 0.088 
1.405 ± 0.159 ±0.168 
1.438 ± 0.206 ± 0.408 
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FIG. 15: Statistical (dashed line), systematic (doted line), 
and total (full line) uncertainties on A^i as a function of the 
parameter a of Eq. (|64p . 



a linear combination of the uncorrected asymmetries 



A' = A- 



(64) 



and choosing the coefficient a in order to minimize the 
total uncertainty on the vahie of A^j . 

It is shown in Sees. IVT -I XIII that the contributions from 
background sources in Eqs. ([9]) and (fT2|) are of the same 
order of magnitude. On the other hand, the dependence 
of A and a on the asymmetry A^^, according to Eqs. ([55)) 
and ([58[) . is significantly different, with Cb ^ Cb- As a 
result, we can expect a reduction of background uncer- 
tainties in (l64l) for a « 1 with a fimited reduction of the 



statistical sensitivity on ^^j. Figure (TSj shows the sta- 
tistical, systematic, and total uncertainties on A^j as a 
function of the parameter a. The total uncertainty on 
Agj has a minimum for a — 0.959, and the corresponding 
value of the asymmetry A'^^ is 

v4^i = -0.00957 ± 0.00251 (stat) ± 0.00146 (syst). (65) 

This value is our final result for A^^. It differs by 3.2 
standard deviations from the standard model prediction 
of Agj given in Eq. ^ . The different contributions to the 



total uncertainty of A^j in Eq. 



65|) are listed in TablelHIl 



The contributions to the uncertainty on A^j for this mea- 
surement are also listed in Table IXIII 



The results ([62[) and (1631) represent two different mea- 
surements of Agj. The uncertainties in Eq. ([5^ are much 
larger because the asymmetry a^ is divided by the small 
coefficient Cb- Since the same background processes con- 
tribute to the uncorrected asymmetries a and A, their 
uncertainties in Eqs. (j62p and (j63p are strongly corre- 
lated. We take advantage of this correlation to obtain a 
single optimized value of A^^ , with higher precision, using 



XV. CONSISTENCY CHECKS 

To check the stability of the result, we repeat this mea- 
surement with modified selections, or with subsets of the 
available data sample. Changes are implemented in a 
variety of tests: 

• Test A: Using only the part of the data sample cor- 
responding to the first 2.8 fb^^. 

• Test B: In addition to the reference selections, re- 
quiring at least three hits in muon wire chamber 
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layers B or C, and the x^ for a fit to a track seg- 
ment reconstructed in the muon detector to be less 
than 8. 

• Test C: Since the background muons are produced 
by decays of kaons and pious, their track param- 
eters measured by the central tracker and by the 
muon system are different. Therefore, the fraction 
of background strongly depends on the x^ of the 
difference between these two measurements. The 
requirement on this x^ is changed from 40 to 4 in 
this study. 

• Test D: The maximum value of the transverse im- 
pact parameter is changed from 0.3 to 0.05 cm, and 
the requirement on the longitudinal distance be- 
tween the point of closest approach to the beam 
and the associated interaction vertex is changed 
from 0.5 to 0.05 cm. This test serves also as 
a cross-check against the possible contamination 
from muons from cosmic rays in the selected sam- 
ple. 

• Test E: Using only low-luminosity events with fewer 
than three interaction vertices. 

• Test F: Using only events corresponding to two of 
the four possible configurations of the magnets, for 
which the solenoid and toroid polarities are identi- 
cal. 

• Test G: Changing the requirement on the invariant 
mass of the two muons from 2.8 GeV to 12 GeV. 

• Test H: Using the same muonp^ requirement, pr > 
4.2 GeV, over the full detector acceptance. 

• Test I: Requiring the muon pr to be < 7.0 GeV. 

• Test J: Requiring the azimuthal angle (/) of the muon 
track to be in the range < < 4 or 5.7 < (/> < 
27r. This selection excludes muons directed to the 
region of poor muon identification efficiency in the 
support structure of the detector. 

• Test K: Requiring the muon ry to be in the range 
\r]\ < 1.6 (this test serves also as a cross-check 
against the possible contamination from muons as- 
sociated with the beam halo). 

• Test L: Requiring the muon rj to be in the range 
|?7| < 1.2 or 1.6 < \r]\ < 2.2. 

• Test M: Requiring the muon rj to be in the range 
|77| < 0.7 or 1.2 < |77| < 2.2. 

• Test N: Requiring the muon r/ to be in the range 
0.7 < |77| <2.2. 

• Test 0: Using like-sign dimuon events passing at 
least one single muon trigger, while ignoring the 
requirement of a dimuon trigger for these events. 



• Test P: Using like-sign dimuon events passing both 
single muon and dimuon triggers. 



A summary of the results from these studies is pre- 
sented in Tables IXIVI and FXVl The last line, denoted as 
"significance" , gives the difference between the reference 
result (column Ref) and each modification, divided by its 
uncertainty, and taking into account the overlap between 
the samples. The statistical uncertainties are used in the 
calculation of the significance of the difference between 
two results. These tests demonstrate an impressive sta- 
bility of the Ag[ result, and provide a strong confirmation 
of the validity of the method. As a result of the variations 
of the selection criteria, all input quantities are changed 
over a wide range, while the asymmetry A^j remains well 
within the assigned uncertainties. For example, the un- 
corrected asymmetry A changes by a factor w 1.5 in test 
C, while the asymmetry A^^ changes by less than 7%. It 
should also be noted that reducing the kaon background 
in test C yields a negative asymmetry A. 

Figure [12] shows the observed and expected uncor- 
rected like-sign dimuon charge asymmetry as a function 
of the dimuon invariant mass. The expected asymme- 
try is computed using Eq. ([T2|) and all the measure- 
ments of the sample composition and of the asymme- 
tries. We compare the expected uncorrected asymmetry 
using two different assumptions for A^j. In Fig. [THla) 
the observed asymmetry is compared to the expectation 
for Agj = 0, while Fig. [TBlb) shows the expected asym- 
metry for Agj — —0.00957. A possible systematic dis- 
crepancy between the observed and expected asymme- 



tries can be observed for A, 



0, while it essentially 



disappears for the measured A^^ value corresponding to 
Eq. (p5|) . It can also be seen that the observed asymme- 
try changes as a function of the dimuon invariant mass, 
and that the expected asymmetry reproduces this effect 
when Agj = —0.00957. This dependence of the asymme- 
try on the invariant mass of the muon pair is a complex 
function of the production mechanism, of the mass of 
the particles being produced and of their decays. The 
agreement between the observed and expected asymme- 
tries indicates that the physics leading to the observed 
asymmetry is well described by the contributions from 
the backgrounds and from decaying b hadrons. 

We conclude that our method of analysis provides a 
consistent description of the dimuon charge asymmetry 
for a wide range of input parameters, even for signifi- 
cantly modified selection criteria. 

In addition to the described consistency checks, we per- 
form other studies to verify the validity of the analysis 
method. These tests are described in Appendices O |D1 
and [E] We determine the asymmetry of track recon- 
struction and the asymmetry of trigger selection. We 
also measure the ratio Fk/ /k using an alternative fit- 
ting procedure. These studies do not show any bias in 
the extracted value of AK . 
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TABLE XIV: Measured asymmetry A^i with reference selections (column Ref) and variations A - H. 



Ref 



B 



C 



D 



ax 10^ 
Ax 10^ 

Q 

[(2 - i^bkg)A - afsS] X 10^ 

Jbkg 



3.731 


1.809 


2.733 


1.809 


1.785 


2.121 


1.932 


1.736 


1.783 


+0.955 


+0.988 


+0.791 


+0.336 


+1.057 


+0.950 


+1.029 


+0.955 


+1.032 


+0.564 


+0.531 


+0.276 


-0.229 


+0.845 


+0.543 


+0.581 


+0.821 


+0.632 


0.959 


0.901 


0.942 


1.089 


1.083 


0.902 


0.915 


1.029 


0.877 


-0.065 


-0.072 


-0.143 


-0.200 


-0.074 


-0.075 


-0.069 


-0.023 


-0.061 


0.409 


0.372 


0.401 


0.303 


0.384 


0.385 


0.426 


0.449 


0.343 



A°i X 10" 

o-(yl*,) X 10^ (stat) 

Significance 



-0.957 
0.251 



-0.976 
0.330 
0.090 



-1.084 
0.293 
0.846 



-0.892 
0.315 
0.324 



-1.107 
0.402 
0.478 



-0.888 
0.328 
0.326 



-1.096 
0.375 
0.498 



-0.873 
0.388 
0.281 



-0.769 
0.336 
0.779 



TABLE XV: Measured asymmetry A^i with reference selections (column Ref) and variations I - P. 



Ref 



I 



K 



M 



N 







N(p^) X 10"" 3.731 2.569 2.208 1.884 

a X 10^ +0.955 +0.896 +1.002 +0.984 

A X 10^ +0.564 +0.407 +0.648 +0.576 

Q 0.959 0.975 0.913 0.895 

[(2 - Fbkg)A - a/s(5] X 10^ -0.065 -0.101 -0.079 -0.125 

Fbkg 0.409 0.439 0.412 0.363 



1.909 
-1.098 
-0.630 

0.877 
-0.142 

0.365 



2.534 
+0.679 
+0.353 

0.940 
-0.081 

0.412 



2.122 
+1.097 
+0.748 

0.949 
-0.019 

0.452 



2.002 
+0.968 
+0.722 

0.983 
-0.046 

0.419 



1.772 
-0.968 
-0.692 

0.934 
-0.044 

0.398 



A^i X 10" 


-0.957 


-1.295 


-0.710 


-0.851 


-0.801 


-0.759 


-1.102 


-0.897 


-0.833 


a{Ali) X 10^ (stat) 


0.251 


0.314 


0.320 


0.320 


0.383 


0.275 


0.344 


0.346 


0.349 


Significance 




1.798 


1.241 


0.482 


0.539 


1.317 


0.622 


0.240 


0.485 



XVI. COMPARISON WITH EXISTING 
MEASUREMENTS 



The measured value of A^^ places a constraint on the 
charge asymmetries of semileptonic decays of B^ and i3° 
mesons, and the CP- violating phases of the i?J; and _B° 
mass mixing matrices. Calculating the coefhcients in 
Eq. ()A9p assuming the current PDG values [2] for all 
parameters (details are given in Appendix E| . we obtain 

All = (0-506 ± 0.043)afi + (0.494 ± 0.043Xi. (66) 

Figure [T7] presents this measurement in the a'^i-a^i plane, 
together with the existing direct measurements of a^j 
from the B-Factories [23 
surement of afi in B2 



Eqs. 
of a'' 



and of our independent mea- 
-J> DsfiX decays [131. Using 
651) and ([66|) and the current experimental value 
0.0047+0.0046 [Hi, we obtain 



Osl 



-0.0146 ±0.0075. 



(67) 



si 



This agrees with our direct measurement of a| 
-0.0017 + 0.0091 [24]. 

An independent method for measuring (ps is through 
B1 — > J/tp4> decays. Such measurements have been per- 
formed by the DO [H and CDF [2^ Cohaborations. All 
measurements are consistent and the combined value of 
(j)s differs from the standard model prediction by about 
two standard deviations [27| . 

Taking into account the experimental constraints on 
Fig. HH shows the 68% and 95% C.L. regions 
obtained from our measurement. The 



68% and 95% C.L. regions from the DO measurement 
using the i3° — > J/V'0 decay [25] are also included in this 
figure. Since the sign of AVg is not known, there is also 
a mirror solution with (pg — > — tt — (f>s, corresponding to 
the change AFs — > — AFg. It can be seen that the DO 
results are consistent. Figure [12] shows the probability 
contours in the ((^igjAFs) plane for the combination of 
our measurement with the result of Ref. uM . 



a, 



si 

of AF, and 



XVII. CONCLUSIONS 



We have measured the like-sign dimuon charge asym- 
metry Agj of semileptonic 6-hadron decays: 



A^i = -0.00957+0.00251 (stat) ± 0.00146 (syst). (68) 



This measurement is obtained from a data set corre- 
sponding to 6.1 fb~^ of integrated luminosity collected 
with the DO detector at Fermilab Tevatron collider. It is 
consistent with our previous measurement IJjj obtained 
with 1 fb^^ and supersedes it. This asymmetry is in 
disagreement with the prediction of the standard model 
by 3.2 standard deviations. This is the first evidence 
for anomalous CP-violation in the mixing of neutral B- 
mesons. 
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FIG. 18: (Color online) The 68% and 95% C.L. regions of 
probability for AFs and (j)a values obtained from this mea- 
surement, considering the experimental constraints on afi [23 ]. 
The sohd and dashed curves show respectively the 68% and 
95% C.L. contours from the B^ — ^ J/'4'4> measurement [251 ]. 
Also shown is the standard model (SM) prediction for (ps and 

Ar,. 



FIG. 16: The observed and expected hke-sign dimuon charge 
asymmetries in bins of dimuon invariant mass. The expected 
asymmetry is shown for (a) A^i — 0.0 and (b) A^i — —0.00957. 
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FIG. 17: (Color online) Comparison of A'^i in data with the 
standard model prediction for aj and aj. Also shown are 
the existing measurements of afi [23l and afi [2J|. The error 
bands represent the ±1 standard deviation uncertainties on 
each individual measurement. 






0.0 



-0.2 



-0.4 



D0, 2.8 - 6.1 fb'^ 





SM 



68% CL 
95% CL 
99% CL 



0s [rad] 



FIG. 19: (Color online) Probability contours in the ((pa,ATs) 
plane for the combination of this measurement with the result 
of Ref. [23I, using the experimental constraints on afi [23 ]. 
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Appendix A: Theory 

This Appendix is included for completeness and to de- 
fine the notations. Assuming CPT symmetry, the mixing 
and decay of the B^, B^ pair {q — s, d) is described [2^ 
by 



.d_ 



B°{t) 



Mg 

(Mr 

Blit) 
Blit) 



Ma 



r r 

J- q J- q 



(Al) 



where Mq, M^^, Tq, and Fj^ are the elements of the 
mass matrix of the B^Bq system. The matrix element 
Mq^ is due to box diagrams [2| . New particles foreseen in 
extensions of the standard model can contribute to these 
box diagrams, and physics beyond the standard model 
can therefore modify the phase and amplitude of M^"^ . 
The eigenvalues of the mass matrix in Eq. (|A1I) are 



Mq + ^AMq~'-{Tq-^ATq), 



M„ 



1 



1 



AAf,--(F,-|--AF,) 



(A2) 
(A3) 



where, by definition, AMq > 0. Notice the sign conven- 
tions for AMq and AF^. With this convention, AF^ is 
positive in the standard model. A violation of the CP 
symmetry is caused by a non-zero value of the phase 



arg 



Ml_ 



(A4) 



The observable quantities are Mq, Fg, AMq, ATq and 
n, with 



AMq ^2\mI'^\ 



Ar„ = 2 Fi 



C0S( 



(A5) 



The charge asymmetry a^j for "wrong-charge" 
semileptonic i?°-meson decay induced by oscillations is 
defined as 



T{B^q{t)^^^+X)-T{Bl{t)^^,-X) 

T{mt) -> ^l+x) + r{B"{t) ^ fi-x) ' 



(A6) 



This quantity is independent of the lifetime t, and can 
be expressed as 



I q I 

\ q I 






(A7) 



The like-sign dimuon charge asymmetry A^j for 
semileptonic decays of b hadrons produced in proton- 
antiproton (pp) collisions is defined as 






Nr 



N^ 



Nr 



(A8) 



where Nj^~^ and Nj^~ are the numbers of events contain- 
ing two b hadrons that decay semileptonically, producing 
two positive or two negative muons, respectively, with 
only the direct semileptonic decays b — >■ /iX considered 
in the definition of N^^^ and Nr~ . The asymmetry A'^^ 



can be expressed [Ifll as 



where 






Zr, 



Vq 



fdZd + fsZs 



AF„ 



2F/ 

AMa 



F„ 



(A9) 

(AlO) 
(All) 
(A12) 



with q = d, s. The quantities fd and fs are the pro- 
duction fractions for b ^ B^ and b ^ B^ respectively. 
These fractions have been measured for pp collisions at 
the Tevatron 2j: 



fd = 0.323 ±0.037, 

/, = 0.118 ±0.015. (A13) 

All other parameters in (jA9p are also taken from Ref . [2J : 

Xd = 0.774 ±0.008, 

Vd = 0, 

Xs = 26.2 ±0.5, 

Vs = 0.046 ±0.027. (A14) 

Substituting these values in Eq. (jA9l) . we obtain 

^^1 = (0.506 ± 0.043)afi + (0.494 ± 0.043Xi. (A15) 

Using the values of af^, al^ from Ref. jlj], 

af,(SM) = (-4.8l}:°) X 10-4 
<i(SM) = (2.1±0.6) X 10-^ (A16) 

the predicted value of A\^ in the standard model is 



4VSM) = (-2.3 



+0 j) X 10- 



(A17) 



The current experimental values of the two semileptonic 
asymmetries are ai_= -0.0047 ± 0.0046 [23|] and al^ = 
-0.0017 ±0.0091 Q. 

It can be concluded from Eq. (|A17|) that the standard 
model predicts a small negative value of v4gj with rather 
small uncertainty. Any significant deviation of A^j from 
the SM prediction on a scale larger than that of the un- 
certainty on ^g[ , would be an unambiguous signal of new 
physics. 

The asymmetry A^^ is also equivalent to the charge 
asymmetry of semileptonic decays of b hadrons to "wrong 
charge" muons that are induced by oscillations [10[, i.e.. 



b ^ riB^^i+X)-T{B^fi-X) ^ , 
"'''" F(B^/i+X)+F(B^^-X) '''■ 



(A18) 
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Without initial flavor tagging it is impossible to correctly 
select the decays producing a muon of "wrong charge" 
from a sample of semileptonic decays of b quarks. The 
charge asymmetry of semileptonic b hadron decays, con- 
trary to the charge asymmetry of like-sign dimuons, is 
therefore reduced by the contribution of decays produc- 
ing a muon with "correct" charge, and is consequently 
less sensitive to the asymmetry A'^^. 

New physical phenomena can change the phase and 
magnitude of the standard model Ml"^'^^ to 



M. 



where 



12 _ 



^,^12,SM . ^^ ^ ^,jl2,SM . I^^l ^^4>f ^ (^19) 



iSM 



= 0.0042 ±0.0014. 



(A20) 



Other changes expected as a result of new sources of CP 
violation jl| are 

AAf, = AA/f^^-|A,| = (19.30±6.74) ps^^-jA^I , (A21) 



Ar, =2\Tf\cos(l)s = (0.096 ±0.039) ps^^ -cos^^, 



Ar, 

AM, 




(4.97±0.94)-10- 



(4.97 ±0.94) -10- 



(A22) 

cos 0s 

(A23) 
sin (A, 



(A24) 

The i?" — > J/ipcj) decay can also be used to investigate 
CP violation. In that case, the CP-violating phase ob- 
tained from fits to the B^ — >■ J/ip'P angular distributions 
is modified as follows [l|: 



0W = -2/3f'^ + 0f, (A25) 

where 

/3f ^^ = aTg[-Vt,Vtl/{VasV,l)] = 0.019 ± 0.001 (A26) 

and the quantities Vts, Vtb, Vcs, and Vcb are the parame- 
ters of the CKM matrix. The contribution of new physics 

to 0s and (j)s are identical. 

Appendix B: Reconstruction of exclusive decays 
1. Reconstruction of Ks mesons 

The Ks meson is used to reconstruct the K*^ me- 
son |11| and to measure the fraction and asymmetry of 
TT — )■ /Lt tracks. The Ks — >■ tt+tt" decay is reconstructed 



by requiring two tracks with opposite charge. Each track 
must have an impact parameter significance with respect 
to the interaction vertex > 3, where the significance is de- 
fined as y^[eT/cr{eTW + [^L/cr{<^L)Y, and er (ei) is the 
projection of the track impact parameter on the plane 
transverse to the beam direction (along the beam di- 
rection), and a{eT)[(j{€L)\ is its uncertainty. At least 
one of the tracks must have an impact parameter signif- 
icance > 4, and at least one of the particles must have 
Pt > 1.5 GeV. The two tracks must share a common ver- 
tex that is separated from the primary interaction point 
by more than 4 mm in the transverse plane. The signifi- 
cance of the reconstructed impact parameter for the Ks 
must be < 4. All Ks candidates satisfying these selec- 
tion criteria are used to reconstruct the K*^ — >■ Ks'rr'^ 
decay. In addition, for the measurement of the fraction 
and asymmetry of tt — > /i tracks, we require that one of 
the pions from Ks decay pass the muon selection given 
in Sec. HTH 

Figure [201 displays the tt+tt^ invariant mass distribu- 
tion of Ks —7- TT^^TT^ candidates in the inclusive muon 
sample for all tt — > /i with 7.0 < px < 10.0 GeV. 
We show separately the sum and the difference of the 
distributions for the samples with positive and negative 
TT ^ /i tracks, which are used to measure the asymmetry 
of TT ^- /i tracks. The Ks signal is fitted with a dou- 
ble Gaussian, and the background is parameterized by 
a third degree polynomial for the sum of the two distri- 
butions and a straight line for their difference. While 
fitting the difference of the distributions all the param- 
eters describing the Ks signal, except its normalization, 
are fixed to the values obtained from the fit to the sum 
of the distributions. 



2. Reconstruction of K* mesons 

The if *+ [ll| signal is obtained by combining the re- 
constructed Ks meson with an additional track which 
is assigned the mass of the charged pion. The Ks can- 
didate must satisfy the track selection criteria given in 
Sec. IIIH except for the requirements on the number of 
hits in the tracking detectors and the x^ of the track 
fit. The invariant mass of the tt+tt^ system must be 
480 < A/(7r+7r") < 515 MeV. The additional track must 
have at least 2 axial and 1 stereo hits in the silicon mi- 
crostrip detector, at least 3 axial and 3 stereo hits in the 
fiber tracker, and a track impact parameter significance 
< 3 relative to the interaction vertex. The cosine of the 
angle between the direction of the Ks meson and the 
additional track must be greater than 0.3. The Ks and 
the additional track must be consistent with sharing the 
same interaction vertex. 

Figure [21] shows the Kstt^ invariant mass distribution. 
The K*^ signal is fitted with a relativistic Breit-Wigner 
function convoluted with a Gaussian resolution, and the 
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FIG. 20: The ■k'^-k~ invariant mass distribution of Ks can- 
didates in the inclusive muon sample for tt — >■ /i with 
7.0 < pT < 10.0 GeV. (a) The sum of the distributions for 
positive and negative vr — > /i tracks and (b) their difference. 
The solid lines present the result of the fit; the dashed lines 
show the background contribution. 



FIG. 21: (a) The Kstv^ invariant mass distribution of K*'^ 
candidates in the inclusive muon sample. The solid line 
presents the result of the fit; the dashed line shows the back- 
ground contribution, (b) The difference between data and the 
fit result. 



background is parameterized by the function 

/bck(M) = {M-Mk-M^Y° 

X exp(piM + paM^ + p^M^ 



(Bl) 



which includes a threshold factor. Here M is the Ks'rr'^ 
invariant mass, and poi Pij P2 and p^ are free parameters. 
Figure I^TTb). which shows the difference between data 
points and the result of the fit, demonstrates the good 
quality of the fit with x^/d.o.f. = 54/49. The measured 
width of the K*+ meson is T{K*+) = 47.9 ± 1.4(stat) 
MeV, which is consistent with the current PDG value [2|- 



3. Reconstruction of K* mesons 

The K*^ meson is reconstructed by selecting two tracks 
of opposite charge and assigning one of them the mass of 
the charged kaon. This particle is required to be iden- 
tified as a muon and to pass the muon selection criteria 
given in Sec. IIIII The second track is assigned the mass 
of a pion and required to satisfy the criteria used to select 
the pion in the K reconstruction. 

Figure [22] shows the K^it~ [11| invariant mass distri- 
bution of the K*^ candidates with K ^f ji in the inclu- 
sive muon sample, while Fig. 1231 shows the corresponding 
mass distribution in the like-sign dimuon sample. 



The measurement of the number of K*'^ — > A'+tt^ de- 
cays with K -^ fi is complicated because of the large 
combinatorial background under the K*'^ signal, and be- 
cause of the contribution of light meson resonances de- 
caying to 7r+7r~. The most important contribution comes 
from the p*^ — > tt+tt" decay with tt — > /i. It produces a 
peak in the mass region close to the K*'^ mass. Figure [24l 
shows the mass distribution of simulated p*^ — )■ tt+tt" de- 
cays with one pion assigned the kaon mass. This pion is 
also required to satisfy the track selections. 

To overcome these misidentification difficulties, the fit 
to the K'^TT~ mass distribution is performed in several 
steps, assuming for the width of the K*'^ meson the value 
obtained in the previous section for the K*~^ meson. Con- 
trary to the K~^7T~ system, the decays of light resonances 
do not contribute into the KsTr~^ mass distribution be- 
cause the Ks meson is identified unambiguously, and the 
K*'^ signal is clean and unbiased. 

The K*'^ mass and detector resolution for K*'^ — > 
AT+TT" with K —>■ n are obtained from a fit to the differ- 
ence of the K^TT~ and A'^7r+ mass distributions. The 
reconstruction efficiency oi K ^ fi tracks demonstrates 
a charge asymmetry w6%, (see Sec.lXI|l. Such an asym- 
metry is significantly smaller for tt — > /i tracks, and the 
contribution of p — >■ tt+tt" and other light resonances is 
therefore suppressed in the difference of the AT+tt" and 
K^TT^ mass distributions. In addition, the contribution 
of the combinatorial background is significantly reduced, 
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FIG. 22: (a) The K'^tt~ invariant mass distribution of i^*" 
candidates in the inclusive muon sample. The solid line cor- 
responds to the result of the fit and the dashed line shows the 
contribution from the combinatorial background. The shaded 
histogram is the contribution of p" — > -k'^-k~ events, (b) The 
difference between data and the result of the fit. 



FIG. 23: (a) The K'^iv' invariant mass distribution of K*'^ 
candidates in the like-sign dimuon sample. The solid line 
corresponds to the result of the fit and the dashed line shows 
the contribution from the combinatorial background. The 
shaded histogram is the contribution of p" — >■ tt^tt" events, 
(b) The difference between data and the result of the fit. 



as can be seen in Fig. ISSl a). The K*'^ signal is fitted 
with a relativistic Breit-Wigner function convoluted with 
a Gaussian resolution, and the background is parameter- 
ized by the function (jBip . Figure ISST b). which shows the 
difference between data and the result of the fit, indicates 
a moderate quality for the fit, with x^/d.o.f. = 71/52. 
The fit gives a{M) = 12.2 ± 1.5(stat) MeV for the K*^ 
mass resolution of the detector. The mass difference be- 
tween K*^ and K*+ is 



M{K*'') ~ M{K*+) = 3.50 ± 0.66(stat) MeV, (B2) 
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which Js consistent with the PDG value of 4.34 ± 0.36 

MeV pt 

The number of K*'^ events in the inclusive muon and 
in the like-sign dimuon samples is determined from the 
fit of the mass distributions shown in Figs. [52] and [531 
The signal is parameterized with the convolution of a 
relativistic Breit-Wigner function and a Gaussian resolu- 
tion. The A'*" mass and width and the detector resolu- 
tion are fixed in the fit to the values obtained from the 
fit to the distribution in Fig. [25k . The mass distribution 
of the p^ — ?► 7r+7r~ background is taken from the MC 
simulation. To improve the quality of the fit, the param- 
eterization of the background is modified by adding two 



FIG. 24: The "_R'"'"7r~" invariant mass distribution of simu- 
lated p'^ — >■ ■k'^-k" decays, where the mass of the charged kaon 
is assigned to one of the two reconstructed tracks. 



additional Gaussian terms: 



P4 exp(- 



(M - Mif 
2al 



+ P5exp( —^ ). 



2al 



(B3) 
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FIG. 25: (a) The difference of the K'^-k' and A"~7r+ mass 
distributions of K'^ candidates in the inclusive muon sample. 
The solid line represents the result of the fit, while the dashed 
line shows the background contribution, (b) The difference 
between data and the result of the fit. 



Here fhckiM) is given in (|Bip . The additional terms are 
needed to describe the distortion of the smooth behav- 
ior of the combinatorial background at large masses of 
M sa 1.15 GeV, due to the contributions of other light 
resonances, and should be considered as a parameteriza- 
tion of the observed mass distribution rather than the 
contribution from specific sources. These terms are sig- 
nificant only because of the large statistics of the inclu- 
sive muon sample, vi^hich contains about 10^ entries per 
bin in Fig. [521 The fit of the K*° signal in the like-sign 
dimuon sample is of the same quality without these ad- 
ditional terms. The results do not change significantly 
if these terms are omitted. The impact of these terms 
on the final measurement is included in the systematic 
uncertainties of the fractions fx and Fk discussed in 

Sec. lyml 

In the fit to the inclusive muon distribution, the pa- 
rameters Ml, M2, CTi, and cr2 and the contribution of 
/5° — >■ TT+TT", are allowed to vary. The ratio of the frac- 
tions of p° and K*'^ mesons is constrained within 10% 
of the value obtained in the simulation. The fit yields 
Ml = 1.095 ± 0.005 GeV and M2 = f .170 ± 0.007 GeV, 
which is far from the region of the K*'^ mass, and does not 
influence the fitted number oi K*° mesons. The x^/d.o.f. 
of the fit is 81/45. The results of the fit and correspond- 
ing residuals are shown in Fig. [211 

In the fit to the like-sign dimuon distribution the pa- 
rameters Ml, M2, (Ti, and 172 are fixed to the values 



obtained in the fit of the inclusive muon sample. The 
xVd-O.f. of the fit is 48/52. The results of the fit and 
corresponding residuals are shown in Fig. [23| 



4. Reconstruction of 0(1020) mesons 

The (/)(1020) meson is reconstructed by selecting two 
tracks with opposite charge, and assigning both of them 
the mass of the charged kaon. One track is required to 
pass the track selections of Sec. IIIII The second one sat- 
isfies the same selection criteria as the pion in the K*^ 
and K*'^ reconstructions. 

Figure [551 shows the K^K^ invariant mass distribu- 
tion of the 0(1020) -^ K^K~ candidates in the inclu- 
sive muon sample, with an additional requirement on 
the transverse momentum of the kaon misidentified as 
a muon, 4.2 < pT < 7.0 GeV. We display separately 
the sum and the difference of the distributions for the 
K -^ fi tracks with positive and negative charges, as 
done in the case of Ks candidates, and use these distri- 
butions to measure the asymmetry for K ^- fj, tracks. 
The 0(1020) signal is fitted with a double Gaussian and 
the background is parameterized by the threshold func- 
tion 



/bck(M) = {M-2MKr" 

X exp{piM + P2M'^ 



'Pz^r 



(B4) 



All the parameters describing the signal, except its nor- 
malization, are fixed in the fit of the difference of the 
invariant mass distributions to the values obtained from 
the fit to the sum of the distributions. 



5. Reconstruction of A baryons 



The selection of A — > pn^ decays [ll| follows that of 
Ks —7- TT+TT", except that one of the tracks is assigned 
the mass of the proton. Figure [ITj shows the p7r~ in- 
variant mass distribution of A -> pi:^ candidates in the 
inclusive muon sample, with an additional requirement 
on the transverse momentum of the proton misidentified 
as a muon, for 4.2 < px < 7.0 GeV. Also in this case we 
display separately the distributions for the sum and the 
difference of the distributions for A and A decays, and 
use them to determine the asymmetry for p — >■ /i tracks. 
The A baryon signal is fitted with a Gaussian, while the 
background is parameterized by a fourth (second) degree 
polynomial for the sum (difference) of the invariant mass 
distributions. All parameters describing the signal, ex- 
cept its normalization, are fixed in the fit of the difference 
of the invariant mass distributions to the values obtained 
from the fit to the sum of the distributions. 



28 



> 



.a 

& 



> 



xlO 
6000 

4000 

2000 



D0,6.1fb^ 

%^/dof= 64/27 




98 



15000 






& 



10000 



5000 



D0,6.1fb^ 



1.02 1.04 1.06 

M(K*K") [GeV] 




(b) 



,V***+A 



1.02 1.04 1.06 

M(K*K") [GeV] 



> 

<u 



M 






20000 



S 10000 



600 




D0, 6.1 fb^ 



%7dof= 33/19 



08 1.1 1.12 1.14 1.16 1.18 

M(pV) [GeV] 




1.1 1.12 1.14 1.16 1.18 

M(pV) [GeV] 



FIG. 26: The K+K' invariant mass distribution of 0(1020) 
candidates in the inclusive muon sample for K ^^ fi with 
4.2 < pt < 7.0 GeV. (a) The sum of distributions of positive 
and negative K ^ ^, tracks and (b) their difference. The 
solid lines present the result of the fit and the dashed hnes 
show the background contribution. 



FIG. 27: ThejJTT" invariant mass distribution of A candidates 
in the inclusive muon sample for p — > /i with 4.2 < pr < 7.0 
GeV. (a) The sum of distributions of positive and negative 
p -^ fi tracks and (b) their difference. The solid lines present 
the result of the fit and the dashed lines show the background 
contribution. 



Appendix C: Track reconstruction asymmetry 



The asymmetry a is defined in Eq. ([6]), and the asymme- 
tries Oopp and Coquai are computed as 



In this measurement of A^j, we assume that the charge 
asymmetry of track reconstruction cancels as a result of 
the regular reversal of the magnets polarity. The method 
developed in Sees . FVl and IXII provides a way of evaluating 
this assumption in data by comparing the track charge 
asymmetry attack with that expected from already anal- 
ysed sources. We select events with one reconstructed 
muon and at least one additional track satisfying the se- 
lection criteria of Sec. IIIII The muon is not used in this 
study, but is required since the events are collected with 
single muon triggers. Nevertheless, the charges of the 
muon and the additional track can be correlated, and 
the asymmetry of the inclusive muon events can bias the 
observed track asymmetry. To eliminate this bias, we 
consider separately events in which the muon and the 
track have equal and opposite charges. The asymmetries 
in these samples can be expressed as 



a-equal 



^track 
^track 



(CI) 



^opp 



'^equal 






+ n 
— iir 



(C2) 



where n*^ is the number of events containing a track with 
charge i and a muon with charge j . The track asymmetry 
fltrack is computcd as 



^track 



„ (Oopp 



equal j- 



(C3) 



As discussed in Sec. IXH we expect the value of attack 
to contain a contribution from the asymmetry of kaon re- 
construction, even after averaging over the different mag- 
net polarities. The expected value of atrack is therefore 
given by: 



atrack 



track /"track 
''K Jk j 



(C4) 



where f^'^^^ is the fraction of reconstructed tracks that 
are kaons, and a^^'^'^ is the charge asymmetry of kaon 
reconstruction. 

The fraction f^'^'^^ is measured using K*° — )• K~^Tr~ 
decays [1 1| . The selected charged particle is assigned the 
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kaon mass and combined with an additional track to pro- 
duce the K*'^ candidate. The K*'^ selection and fitting 
procedure is described in Appendix |B] The measured 
number of K*'^ mesons is converted into the number of 
kaons using a method similar to that presented in Sec.FVl 
The measured f^^''^ fraction is assigned a systematic un- 
certainty as in Sec. IVIIII 

The same K*'^ — > K^tt^ decay is used to measure the 
kaon reconstruction asymmetry. This measurement can 
be biased by the asymmetry of the muon, because the 
charge of the kaon and the muon can be correlated. The 
kaon asymmetry is therefore measured separately in the 
samples in which the muon and the kaon have equal or 
opposite charges. The asymmetry a^^^"^^ is computed in 
a way similar to Eq. (|C3|) : 



track _ ^ (r,°PP 



cquaK 



(C5) 



The K*^ mass distribution is plotted separately for posi- 
tive and negative kaons in each sample, and the sum and 
the difference of these distributions is fitted to extract 
the quantity A, corresponding to the difference in the 
number of K*^ decays with positive and negative kaons, 
and the quantity E, corresponding to their sum. The 
asymmetry a^'' is measured as 0°^"^ = Aopp/Sopp, and 
a similar relation is used to obtain the asymmetry a'^^^ . 
The expected and observed track reconstruction asym- 
metry, i.e., the right and left side of Eq. (jC4l) . are com- 
pared in Fig. [28] as a function of track pT- There is 
excellent agreement between these two quantities. The 
X^/d.o.f. for their difference is 5.4/5. The fit of (5track — 



^track ^ 



track ^track 
K 



fj^ to a constant yields the following es- 



timate for a residual track asymmetry of 



-'track 



= +0.00011 ±0.00035. 



(C6) 



We conclude that the residual track asymmetry is consis- 
tent with zero as expected. The uncertainty on this value 
is about a factor of 16 times smaller than the observed 
charge asymmetry in the like-sign dimuon events. This 
study provides an additional confirmation of the validity 
of the method used in this analysis. 



Appendix D: Trigger asymmetry 

We determine the kaon, pion and proton charge asym- 
metries using events passing at least one single muon 
trigger, and we apply the same asymmetries to the like- 
sign dimuon events collected with dimuon triggers. Sim- 
ilarly, we measure the muon reconstruction asymmetry 
using events passing the dimuon triggers, and we use the 
same quantity for events collected with the single muon 
triggers. If the trigger selection is charge asymmetric, 
and this asymmetry is different for single muon and for 
dimuon triggers, the obtained value of A\y can be biased. 
Since we extract the asymmetry A\-y from the difference 
A — aa, and because the value of a is very close to unity, 
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FIG. 28: (a) The expected (points with errors) and measured 
(histogram with negligible uncertainties) track reconstruction 
asymmetry, (b) The quantity 5track = atrack - a^'"''^/x''*. 



our measurement is especially sensitive to a difference be- 
tween the charge asymmetry of dimuon and single muon 
triggers. 

To examine the impact of this difference on our re- 
sult, we repeat the measurement of A^j using dimuon 
events passing any single muon trigger without requiring 
dimuon triggers. The result of this test, given in column 
of Table IXV[ does not indicate a bias from the trig- 
ger selection. In addition, we measure A\-y using dimuon 
events passing both single muon and dimuon triggers. 
The result of this test is given in column P of Table IXVI 
These two tests provide a residual difference 5t between 
the asymmetry of dimuon triggers and single muon trig- 
gers of 



5t = +0.00010 ±0.00029. 



(Dl) 



From this result, we conclude that the trigger selections 
do not produce any significant bias to the value of A^j. 



Appendix E: Alternative measurement of Fk/ fx 

We consider two methods to obtain the parameters 
Oi = F^/ f^j^ in the px interval i from 



a. = 2 



Ni{K*'^ -^ K ^ ^i) 



Hi 



N, 



ni{K*" -^ K ^ fi) 



, (El) 



where Ni refers to the like-sign dimuon sample and rii 
refers to the inclusive muon sample. Ni{K*^ —>■ K —i' 
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TABLE XVI: Values of a, = Fk/fk obtained through two 
methods, with their statistical uncertainties. 



Bin 


Q, from Table ITTTl 


Oi from null fit 



1 
2 
3 

4 


1.309 ± 0.340 
0.987 ± 0.082 
1.022 ± 0.050 
0.758 ± 0.101 
1.406 ± 0.159 


0.954 ± 0.217 
0.942 ± 0.069 
1.031 ± 0.027 
0.806 ± 0.055 
1.292 ± 0.079 


All 


0.998 ± 0.038 


0.990 ± 0.022 



/J,) and ni{K*'^ -^ K ^ fj) are obtained by fitting the 
invariant mass histograms of K*^ — > ttK — > ji in the like- 
sign dimuon and inclusive muon samples, respectively. 
These fits require precise modeling of the p" resonance 
and of other backgrounds. In the first method we use 
the results listed in Table IIIII and obtain the values of a; 
listed in Table EVB 

A second set of values for the ai parameters can be 
obtained by finding a scale factor which minimizes the 
differences between invariant mass distributions for the 



K*^ candidates in the inclusive muon and dimuon sam- 
ples (null fit method). Invariant mass distributions for 
A'*" candidates analogous to the ones shown in Fig. [22] 
and[221are built for each bin oi K ^ ^ transverse momen- 
tum. We scale the invariant mass distribution in the in- 
clusive muon sample by a factor aiNi/(2ni) and subtract 
it from the invariant mass distribution obtained from the 
dimuon sample. The contributions from the p resonance 
and from other backgrounds cancel to first order in the 
difference of the two invariant mass distributions, sim- 
plifying the convergence of mass fits. We then vary the 
factor ai and find the value that yields a null normal- 
ization factor for the residual K*'^ signal. The statistical 
uncertainty on a.i is obtained by choosing the value of 
the scale factor which yields a positive (negative) value 
of this normalization factor which is equal to its uncer- 
tainty. The values of ai obtained with this method are 
also listed in Table IXVH with their statistical uncertain- 
ties. There is remarkable agreement between the values 
of ai obtained using the two methods. 
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